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PREFACE 


The wave drag reduction procedure presented in Volume I of this report makes 
use of three major computer programs: the Initial Value Surface (IVS) program, 
the Three-Dimensional Method of Characteristics (3DMoC) program, and the Surface 
Fit and Minimum Search program, which are described in the present volume in 
three parts. The procedure can be used to improve an existing aircraft or to aid in 
the design of a new one, In the case of improvement, the existing aircraft naturally 
serves as the baseline, In the case of new design, the preliminary configuration, 
which is usually obtained through considerations other than the wave drag, can serve 
as the baseline for wave drag reduction. 


The wave drag reduction procedure consists of three major steps: 


1, The first step is to select either the 5 x 5 Latin Square for the forward 
fuselage or the 3 x 3 Latin Square for the blended wing configuration and to 
describe the baseline and variation configurations using the body description 
method presented in Appendix C of Volume I. The sections on the descrip- 
tion and variation of configuration in Appendix A or the main text of Volume 
I should be consulted in producing the body descviption. This step is time- 
consuming but must be done carefully to assure success in wave drag cal- 


culations by the method of characteristics. 


2. For the next step of calculating the wave drag coefficients of the configura- 
tions sampled by each cell of the Latin Square, the IVS and 3DMoC programs 
are used. These two programs can run together with the IVS program provid- 
ing the initial data needed by the 3DMoC program. Input cards can be pre- 
pared according to the instructions given in Parts 1 and 2 of the present volume 
and fed to the IVS and the 3DMoC programs for wave drag calculations. Care 
needs to be taken in the preparation of these cards, for if the calculation fails 


to proceed further the first item to check is the correctness of the input cards. 


3. After the wave drag coefficients have been calculated by the 3DMoC program, 
the Surface Fit and Minimum Search Program can be used to define the mini- 
mum wave drag configuration for a given set of constraints as explained in 
Part 3 of the present volume. The program can be used to calculate various 


minimum wave drag configurations for different sets of constraints as long 


| Weceding Fe9e BLavh - FiLMEL 


as the same baseline configuration is considered and the variables stay 
within the ranges. 


In the following parts, detailed information for setting up and using the three 
computer program are presented. Samples of deck set-up are compiled in the 
Appendix for illustration and reference. 
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PART 1: THE INITIAL VALUE SURFACE PROGRAM 


SECTION I 


INTRODUCTION 


The Initial Value Surface Program is used to generate starting data for the Three- 
Dimensional Method of Characteristics Program. In combination, these two programs 
determine the supersonic flow past smooth blunted bodies for an angle of attack 
at which the initial value surface remains on the blunted nose. Since the sonic point 
on a sphere lies between 40 and 50 degrees from the stagnation point, the sphere-body 
junction should occur at an angle greater than approximately (50 -a@) degrees on the 


sphere, measured in the body axis system. 


The lowest usable Mach number is around 2.5, with lower values attain- 
able by careful attention to the initial value surface solution. There is no computa- 
tional upper limit on the Mach number. However, above Mq=7.0 (approximately) 


the ideal gas model becomes somewhat inaccurate, 


The Initial Value Surface Program consists of four separate programs, all con- 


nected through an OVERLAY sequence. These programs are: 


(1) Blunt Body Program 

(2) Initial Value Line Interpolation Program 

(3) Rotationally Symmetric Method of Characteristics Program 
(4) Initial Value Surface Interpolation Program 


The Blunt Body Program determines the subsonic through slightly supersonic flow 
over the blunted nose. The Initial Value Line Interpolation Program interpolates in 
the Blunt Body Program data to develop a start line for the Rotationally Symmetric 
Method of Characteristics Program. This latter program calculates the required re- 
maining supersonic flow field over the nose. The Initial Value Surface Interpolation 
Program interpolates in the RSMoC-generated data to determine properties on a plane 
normal to the body axis, given the number of data rings and the number of meridian 


planes. 


The final output is provided in printed form and is written on a scratch tape for use 


by the Three-Dimensional Method of Characteristics Program. 


The Initial Value Surface Program consists of approximately 3000 cards. It is 
written entirely in FORTRAN EXTENDED and has been run under OP'T = 2 on the 
CDC FTN Compiler. 


U 


ch stot eh 808 enn tt eri tte te ee gh Rt 


_Fepoant “ , i, i Grn : inte =s 


SECTION II 


INPUT/OUTPUT DESCRIPTION 


1. INPUT DESCRIPTION* 
The input data cards required by the Initial Value Surface Program are shown in 
Figure 1. They are described in detail below. The limits indicated refer to storage 


limitations of the program in its present form. 


Card 
No. Variable Format Description 
@ [. TITLE 7A10 Problem identification 
ACHM E12,4 Freestream Mach number : 
GAMA E12.4 Ratio of specific heats 
ALFA E12.4 Angle of attack, in degrees 
; ELAM E12.4 Slope of the sphere at the sphere-body juncture, 
in degrees. For bodies with non-conical shapes, 
this slope should be the largest value of the 
slopes around the sphere-body juncture. 
RINGS E12.4 Number of data rings on the initial value surface. 
@ A minimum of 3 is recommended. 
Limit: RINGS < 20 
ONE05 E12.4 Mpody> Mach number on the body at the point 
¢ through which the IVL is passed. For Mg22.5, 


a value of 1.05 is recommended for Mpody: 


creased to 1.10. The upper limit is the highest 
Mach number calculated by the Blunt Body 
portion of the program on the surface of the 


x For lower Mach numbers, Mpody should be in- 
| 
| 

sphere. When Mpody is specified to be greater 


‘ than this highest value, the program automati- 
cally reduces Mpody to .01 below this value. 


*The Appendix may be consulted for sample deck setup. 
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Set dk RR 


Card 
No. Variable 


BBW 


BTEST 


STEPS 


POINTS 


EM 


Format 


E12.4 


E12.4 


E12.4 


E12.4 


E12.4 


E12.4 


@ if B=0, omit card 4 


G 


DETA 


© FJ 


E12.4 
E12.4 


E12.4 


Il 


F3.0 
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Description 


Body bluntness, Bhody- For a sphere, Bhody 
is equal to 1.0, For an ellipse, Bhody is equal 
to a?/b?, where a and b are the semi-axes of the 
ellipse. 

The tolerance, AB, accepted in Bpody- A value 
of .001 is usually sufficient. : 

Number of steps required to reach the body, not 
counting the shock wave. For Mw22.5, a value 


of 4 should be sufficient. For lower Mach num- 
bers, larger values should be used. 


Limit: STEPS < 10 


Initial guess at the shock wave bluntness, Bepock: 


If B=0, the program computes its own value 
based on an approximating curve. Card 4 is 
not required for this case. 


Number of shock points. For M@>2.5, a value 
of 70 should be sufficient. For lower Mach 
numbers, this value should be increased. 


Limit: POINTS <100 


Number of meridian planes in 90°. This defines 
the spacing of the meridian planes. Note that, 
if EM is the number of meridian planes in 90°, 
2EM-1 is the number in 180°. 


Limit: EMs 10 


Step size multiplier, to be obtained from previous 
calculations or other knowledge. 


Initial step size normal to the shock, AN 


Initial step size along the shock, A § 


Output control for the RSMoC portion of the 
program. I=2 is recommended. 
=1, only body points are printed 
If I | =2, all points are printed 
=4, detailed print is generated 


Initial C-Summary Number for obtaining local 
checkout print. Should normally be zero. 


ou 


' ania | 


Card 
No. Variable Format Description 
(Q) FK F2.0 Final C-Summary Number for obtaining local 
checkout print. Should normally be zero. 
HOL (3) E18.6 Step size along the body for the RSMoC portion 


of the program. A value of 0.1 is recommended. 
Should difficulties arise in the RSMoC portion of 
the program, a small change in this value will 
rearrange the mesh points and often cure the 
problem. 


Figure 2 shows the input data cards required to generate an initial value surface 
for a blunted 25° cone at 15° angle of attack. 


2. OUTPUT DESCRIPTION 


The first item of output is a listing of the input data, as read. Figure 3 is an ex- 
ample of such output, Whenever unusual messages or stops occur later in the program, 
this listing should be examined for correctness. Also note that an interpolation Mach 
number which is higher than the recommended value (for the freestream Mach number 
specified) has been used. This was done to exercise the automatic resetting feature 


of the program, which was discussed in the preceding section. 


Each iteration of the Blunt Body portion of the program generates output, such as 
shown in Figure 4. At the top are given the freestream conditions, followed by the 
governing quantities for the Blunt Body portion of the program, B, DETA, DY and 
G which have all been identified in the preceding section. Following this is a progress 


report showing the number of points on each shock-like line. Note that the program 
| r will drop a minimum of six points in passing from line to line. More points may be 
big dropped due to increasing roughness in the interpolation process. The results of the 
body fit process are listed, where DELTA is the shock stand-off distance, 4/ Rsho ck? 
B(BODY) is the body bluntness, and R(BODY) is the ratio of Rpody/ Rghock* The 7) 
j and € coordinates and the local Mach number are given for each subsonic body point. 
Near convergence, the final body Mach number should be only slightly less than 1.0. 


When convergence has occurred, the complete blunt body flow field is printed out, 
° with the coordinates now normalized by Ryo dy’ Figure 5 shows a typical first page 
of such an output. At the top are given the freestream conditions, and the geometric 
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eeeee 
e @ MACH NO. = 6.8000 
° © GAMMA = 144000 
. * ANGLE OF ATTACK = 15.000 | 
* U * AODY SLOPE = 25.000 
« © NO. OF RINGS = 5.0000 | 
* © INTERPOLATION MACH NOs = 1.1500 
® *® RODY SLUNTNESS = 1,0000 | 
X.) * BLUNTNESS TOLERANCE = .10000E-02 | 
© Te STEPS = 4.0000 | 
e © BISHOCK) = -0. | 
= * NO. OF POINTS = 70.000 | 
e ® NO. OF MERIDIANS = 4.0000 | 
‘ nl 
FIGURE 3. LISTING OF INPUT DATA 
FREE STREAM CONDITIONS 
»ABO000F +01 
1.400000 
PARAMETERS 4= 24514916 +00 
NETA= +2074 3AE-01 
’ Ove *200000F-01 
G= + 100000E +01 
ETA= 1,0000000 NO OF XI 70 
29792562 NO OF AT 62 
29585124 NO OF KI $5 
2 937TARG NO OF XT 49 
; 29110268 NO OF XK] 43 
' 2B952A10 NO OF KI. 5 = 37 
. Poeoevoscervesoronncecoerecvee 
©1031314E +00 
' BI(BONY)=  .AS19421E +00 
z R«(BODY)= 
ooo 
eecre0ee0@ oe PSPSPS OHSHOOHOHOOH HU GOOOHOOOSEHOPGHHEHOEOHOHERECIO SCORES OOOCOHEHSSOHESOCAOOEHOOOFEH EHO 
NO. FTA xI MACH 
1 +893619E +00 ~100000F -01 «203352€-01 
2 28935666 *00 -300000F-01 +483756E-01 
3 «893461E*00 +500000F-0) - 7A7246E-01 
> « +893302F +00 «700000F -0) «10951 TE +00 
5 eADIO9IE*00 +900000€-01 #1405066 +00 
6 +892825E*00 «1100006 +00 + 171634E +09 
Li +892506E°00 «1300006 +00 +202892E +00 
a +892128E*00 +150000F +00 © 234284E +00 
+ 9 28916946E 200 «1700006 +00 + 265825E +00 
‘ 10 ~891203E°00 «1900006 +00 ©297533E +00 
i +890653E°00 »210000F +00 + 329629E +00 
> 12 ~890041E*00 +230000F +00 © 361538E +06 
\4 13 +BA93I6SE*00 +250000€ +00 «39 3BB7E +00 
& 14 »8AB625E*00 #270000 +00 ©426507E +00 
." 15 «AA7B16E*90 +290000F +00 24594 35E +00 
4 16 ~8A6936E*00 +310000E +00 ©492709E +00 
* 17 «885982E +00 +330000F +00 +526374E +00 
18 +BR“9S0E*00 + 350000E +00 «560483E +00 
t 19 ~AR3B83SE*00 «370000 +00 +595096E +00 
4 20 «892632E +00 +390000E +00 +630282E +00 
21 «8A1337E*00 «4100006 +00 2666126600 
22 «A79942E +00 «4300006 +00 +702726E +00 
23 28 7R440E*00 «450000F +00 + 7402006 +00 
26 2876B21E°00 «4700006 +00 + 77T8697TE*00 
‘ 
FIGURE 4. INTERMEDIATE BLUNT BODY OUTPUT 
. 
. i] 
q j 7 
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COMPLETE FLOW FIELD NORMALIZED BY RBODY 


Me 8.800600 
GAMA= 1.400000 
PR +100000E 01 
R= 2 0O000F +o) 
DELTA= +137283€ +00 


x R MACH NO VELOCITY TOT PRES PRESSURE TEMPERATURE OFNSITY DEL TA THETA 


633329E-06 .129R26€-01 23909266 .T7T6BRECOO §=.S5SO4B4F-02 .90)714F 202 .159990F 202 .S63K0TE*O1 2.5AA52 89,446099 
o5700S2E-03 .3894673E-01 23956237 .179759E +00 = .5514465F-02 .9010266+02 .159875E*02 .S63SRIESO) §=7.69219 88,32308 
eISAITIVE-02 ,669)22F-0) 2404466) .183826E+00 .553374F-02 8996526402 215964666 +92 .563528E 901 12.59567 87.205465 
eILOSI4GE-02 .90B7TOE-01 =. 4183654 = 1A9752E600 =.556279F-02 .897595F +02 .159303E02 .563450E901 17,15868 66.06830 
eSLISOLE-02 .1I6AR42E*00 =. 435 75A3 =. 197357ES00 S601 THF -02 .B894R60E*02 .158847E*N2 .563965E0) 21.3374) 86,9718) 
eTOTHOOE=02 .142B0TE*00 =. #566327 =. 2064416200 .5650RGE-02 .B91GSIESO2 .158279E+G62 .563724E 40} 25,08290 83, 85618 
wLOT256E-O1 = 168772E00 «=. @AODTH2 «=. 216799E00 .571030F-02 .BA737BESO2 .157600E 02 .5630S6E*01 286.3873) 82.74159 
eL4PA9IE-O1 L9GTITE*OO «=. SOTLOLS =. 22R23KE 900 =. 57ROGHE-02 ,8B264BE02 =. 156K) 2E 902 .S628T1E%0) 31.26663 81.6282) 
wDA36846E-01 .220701E*00 «=. 53616146 =, 260579E 200 = .S5RH164F-02 .B772TZEC02 .1S5S9IHES02 .S6245PE 901 33.7625R 80.51622 
e22965)E-9) 8. 24664666F +00 -5671563 42536686 +00 § .595632F-02 .B71261F 02 .156914E 912 = S62417E*01 35.85596 79,40578 
o280R23E-01 .272631F +00 25999329 .267368E 600 = .60SR99F-02 .864629E 002 .1S3ROKESO2 .S621G7TECOL 37.66195 78.29706 
2337232E-01 ©.298596E*00 = 66342815 =. 281565E*00 = .617623F-02 .BSTIBHESOS .1S2601E 902 .S561RGBESO} 39.13708 77.19020 
e39ADLIE-01 3245616200 2670030A .296159E 00 6306676 -02 1669556902 .1ST29SEe4e =. S561S20E°0) 40.3755% 76.08535 
2465R99E-01 «3505266 +00 27070626 =.3L1067E*00 = .645106E-02 .841142E°02 .149893E902 .561160E*01 41.38836 74,98265 
eSBAZIIE“OL =. 3TOG9IESIA §«=—«-. PSLOTE «= 326228E 200 =.662021F-02 .B3I2ITIESO2 .148398E*N2 .S560770E*OL 42.20319 73,88222 
eSISOTTEOL 44024655E 200 =. 784437K =. 341552E 200 .67B8S05F-02 .822657ES02 .146812EeN2 .560347E°01 42,84635 72,78618 
eO99163E-01 2 428420F 00 =. 8 2406R3 =. ISTOLTES00 § .697660F-02 .812619E 02 .145139E9M2 .S59ADIESOL £3.33307 71.68864 
e7TATASIE-01 44543856 +00 sB6554R8K  § .37256R8E 600 3. 7L3960GF-02 .B802077E%02 .143382EON2 .559400F S01 43,62777 70,59568 
eRA2ZLOIE=01 480350E*00 =. 907942AR  . 3ABIO9E +00 =. 741463F-02 .7TILOSOE S02 .1461543E9N2 .SSBATOE*OL 43.97667 69.50540 
ePAIITSE-01 = .506315E*00 = 69509269) § .4037K5E00 =. 766363F-02 .779559E202 = 613962REG2 =. SSBIDIESO] 64.05702 68,4)786 
e10A8754E*00 §=.532280F*00 «=. 994 7B73—-419390E 600 =. 793521F-02 .767624E 202 6137639ESN2 .SSTTIOE SOL 4409762 67,33312 
eTEIBR7E CGO 6 558245E20Q 240995192 4 434956E +00 .A230SKE-02 .755267F+02 .135579E S02 .557N69E*01 44,05610 66,.75122 
e131605E*00 = .SR4210F200 1,0851254 4504646600 .ASSIDIF-02 .7425106°02 .133452Ee02 .556387E 90) £3,.96212 65.17219 
21439166600 .610174E00 161316162 .465893E*00 .A901G4F-02 .729373E202 .131262Fen2 8. 555461E90) 43,474336 64,09605 
2 156828E+00 .636139F 200 1.1790079 .4A1228E+00 .928163F-02 .71SB79F S02 .129013E9N2 .554990E4O1 43.57670 63,02279 
2170353E+00 .662104E*00 =1.2273231 .496452E+00 .969526F-02 .702050F +02 .12670RE*02 .S54071E%01 43.23816 61,95239 
2 1A4L699E 200 =.688069F+00 142765899 .5115566400 .101454E-02 .687907E 602 .124350E 902 55370390) 42.99299 60,88483 
2199278E 600 6714034F200 = 1,326H619 6526521600 .106357E-01 6734726202 .121943E 92 =.552782E401 42.57582 59,82005 
e2th703E600 «6 739999F 600 1,37B1282 .5613464E +00 .)IIO9BE-01 26587666902 4119492E 092 551 205F 20) 42.1107] 58,75798 
22307AGE*09 =. 765964E400 144304632 .556013E00 .117523F-01 .663810F +02 .116999E*F2 .S5O7TIEOL 41.64123 57,69856 
224753R8E 600 =. 791929E400 124839271 65705206600 .123881E-01 .62R626F 202 11446RES02 4 S5&91 THEOL 41.18053 56,646162 
2264977E 600 = .817893E+00 1,5385655 .5AGA59E +00 =6130827E-O1 .613735E 002 .LLI90ZE*O2 .54B012E+O1 40.67237 55,58709 
e2PAI117E*O0 =.843RSBE+00 1.5946601 599024400 .139426F=-01 .S9TASSE%O2 610930GF 202 S66 7B1E201 40.13621 54,53480 
2301976E +00 = .AAIRZ3£200 = 1.651O192 .513009E+00 .146749E-01 .SBI9DHE*02 6106679£992 .545475E 401 39,57715 53, 48458 
e321S7TIE*O0 .BOSTABE*OO 1,7101777 .626R10E+00 .15SB7BE=01 .S566012F 202 1104029402 .S44N91ES01 38.99618 52,43626 
2 3419226600 «6923 753E400 = 1.277019R0 .5460423E+00 .165905F-01 .549987E*02 .101357E*02 45426276001 38,39687 51,38955 
2 IH3069E 600 .9GTTIBESND 1.B8317705 .6538G4E400 = .175937E-01 .533R50E*02 .986669E 201 .541963E201 37477670 50, 3442R 
2 186975E 690. F7IABIESOD «1,R949965)  AATOTIESOO .1A9Q9SE=01 .SITHLIESO2 .95960RE*01 .539406F 20) 37.13693 49,30016 
240772466200 .999848F 200 8 1.9599790 .6R0102E*00 .202520F-0f .5OLRTITE C02 6F224639E 00} «=, SITA4GE*O) 36.4266 68,256846 
24313226200 .1025H1E%OL 7.02664639 = .492934E 000 .217373E=01 .484946F 202 9051296001 .535773E901 35641275 9 47,21403 
e455797E 600 .1OSISBESOL 2.0957224 6 705569F 600 = .233839F=0) 6 46B532F 602 =. ATTTO4F SN) .533779E201 35.17R05 46,17134 
e4RALIBIECOO = LOTTSGF SO} 2.1467570 67180036500 .252136F-01 .452091E*02 .950350E+01 .S31A53E*01 34,42917 45,12836 
2507SO7E*OO = LLIOISIESOL 22601117 27307366400 = .277S1GE=01 .435635E 002 =. AP291OESN1 «529 °RGE 201 33.7166) 44, 08464 
eSIGALOE S00 PP29GTESOL 263159660 6 7427K8E 00 .295269F=01 .4191T9IF S02 .79SKTOE SDL .52695RE01 32.99075 43,63967 
eSOTEITE SIO CPISSSSE SOL 223965171 27541006600 = .320745F=01 .402737ESO2 .THBO4IESNIL 4524363E°01 32.25183 41,99291 
eSYPSTIE C00 CTLANGOR SOL 2.467547K9 .TAST3IF COO 863693696 -01 .386319E 602 27404 70E%01 =. 522581E40) 33.49999 40,96375 


FIGURE 5. FINAL BLUNT BODY OUTPUT 


———ee 


L properties of the flow field, Rpody> Boody? and 4/Rpody: ‘ 


| Be Data along each shock-like coordinate line are printed out, starting with the shock 
itself. The solution is located in a Cartesian coordinate system with the origin located 

? at the stagnation point on the shock (i.e., the intersection of the bow shock and the 
stagnation streamline). X and R are the abscissa and ordinate of the point, respectively, 
MACH NO is the Mach number, VELOCITY is V/Vq , DELTA is the local flow direc- 


tion in degrees, and THETA is the local inclination of the shock wave in degrees. 


Figure 6 shows a typical final page output which represents the body. Note that the 
total pressure values are all identical. This is the distinguishing feature of the body 


point output. Also note that the maximum body surface Mach number is 1.0910140,. 


4 


7 


x 8 MACH NO VELOCITY TOT PRES PRESSURE TEMPERATURE OFNSITY OELTA 


eNISBIIE+00 .LISSTAE=G1 = 60166529 .672H25E-02 =. 36AZBHE-02 = .119105F603 .19H3N2E9N2 612957001 28.61060 
e136I50E +00 .366535E-01 «=. 0450707 .20691KE-OL .36R2AGE-02 1189536903 .196261E*N2 .612401E 90) 47,75170 
21376256600 .STTO9SE-01 «= OTSC9G1 345563801 =. 36AZBGF-O2 .LLB6SIE*O} .194100E*02 .621789E 90} AO.50K4S 
eh 390I7TE*OO = .BOBIALE-O1 = 1059335 =. 4A40SSE-O1 =. 369284E-02 .118199E903 .19IBABESO2 .609422E*01 45.20960 
eLALLGOE*OG .103909F 200 = 13SHZ292 -.622560E-01 .36A284F-02 .LLTSOTE*O3 .193606E 902 .607604E901 43.89142 
eIGIRAIE CON .126965E200 «= EAGLO7TO «= PAIOPRE$OR =. 2GAZ2BHE-02 = BIGB4TECO §««193252E 902 .004437E90) 2.56355 
eL4TL20E 200 .1500Q0E*00 1960900 §=.899697TE-01 =. 1682H4F-02 41159526903 .192828E*N2 .601I24E 901 81.27966 
eVS0902E+00 .17T3Q09E00 «=e 2273011 = 1 D3HIBE*OO =. IHBZ2BGF-02 .114913E903 §=.192333E902 .597469E%01 79.89053 
et55233E90G .EPSPABE*HO «=o 25HUHHFI_ A DIETTIRE SOO §=6 3642846F-02 2333732E003 = -IILIHGESO2 .593N7TECO1 78.54725 
eLOOLISE*00 .228933E*00 = 2690067 § .LBIAL2E*00 .368284F-02 1126126003 .1DLLZ7TESN2 = .588152E901 77219965 
216555466400 .261R39E*00 «= 3201493 21655246600 =. 36A284F-02 .110956E903 19041 7E*02 .582700E901 75.84748 
ePISSIESO0 § .264703F*00 «= 353 927K = 159458E 00 = 3642B4F-02 6109367E903 .189633E9N2 65767266901 74.49028 
eLTALITESOO .287S14E*00 «= BABIN7L3 «I 7IGIGECOU =. 34B2B4F-02 .L0T6GTECO3 .IABTTTESO2 .570235E%01 73.17743 
o1AS253E*00 «6.3107 72E+00 eG1S1159 .LATSLIE*OO .368286F-02 .105801E*03 .187B46E*02 563733890) 71.75613 
0192965E+00 . 332969E +00 24474006 =. 201467E+00 =. 36A2A46F-02 41039316093 .196840EC02 .5557238 00) 70.38142 
2201260E*00 .355H00E*00 =. #800695 =, 215529E 200 =, 3H4286E-02 LOL TSIESO3 ,1ASTSBECO2 ,S4TTIOE*O1 68,99616 
e2LO14RECO0 =. 3781STE*00 «= 5131724 = 229HT0E*OO =. 36A2BHE-02 99534626902 .1AGS9TESN2 .539)96E*01 67.60093 
©21963SE*00 =. 400432 +00 25467665 .2639B2E900 =. 36R2B84E-02 .97Z2129E902 .1AIISTECO2 .530184E90) 66.19414 
e229731E*00 44230186 *00 sSHO4RTH .25BLTHE*00 .36AZAGE-02 .947T799ESO2 .182034E*H2 .520472E 201 44.77386 
224044KF 000 = .445205E+0U 015/022 .272577E*00 .36A28GE-02 .922377E*02 .190625E*02 .S10ASRE COL 63.37779 
e2517TILESO0 =. 46 TUASE +90 20512210 eeRTOY9E COO .369286F-02 .A9SATBECO2 .179LZ2TECO2 .500135E0) 61.89319 
ePOVTAQE*OO =. 4A 9543F +00 6875518 e3OR770E*00 =. 36A2A4F-02 .868I912E602 .17ISISES02 .499NIGE*O) 40.40679 
e2TAG2KE S00 SEIS TOE SOO «=a 7240550 B16427E 200 =, 364284E-02 .BI9KTTESO2 .17TSBG2EON2 .477S1BE+01 58290460 
e2R97TSSE COO §=.533251E*00 «=. 7632804 = I3LAPSE*OO =, 36A2BHE-02 .HODISTESO2 .174OG1E%N2 .465383E 401 57237182 
eVOITSIECOY 6554R69E 690 =. BOIUZ TH «4 3HTOIBE*00 $=, 36AR2B4F-02 .TT9IZIESO2 -1TZIZIEC02 .452458E401 55.80248 
eDEASTOE*OO .ST6307TE*00 «=. B44 SHHO «= BHZTISESOO =, 36A2ZBGE-02 .THTLIOESO2 .17TOOTIES02 .439796E201 S4&.1A922 
2333970600 =.597543F +00 68875710 .I7BBORE*OO .3A8286E-02 .7I3B833EC02 .16TBTIE*N2 .425727E*01 $2.57272 
2 ISOI2Z7E*D0 =.62B55)E°00 49331278 = .395428E400 ,36A2AGE-02 .679153E02 .165499E*02 .410366E01 50.79112 
eVTIZIESOO 2639203F #00 = 6 PALOOTH §.412722E000 =. IHAZBHF-02 .642AOTECO2 .162923E*N2 .396583E901 48.97898 
eVAGIIFE SIN §=.ASYTHIE*00 1403343465 4430A96E 00 .369284F-02 .606679E02 .16009TE*N2 .3ITTAISE*OL 47.04582 
e403611E200 .6798BRE +00 2,09101460 44502436400 ,364284E-02 5641536902 .156955Een2 .359435E+01 45,02390 


FIGURE 6. FLOW PROPERTIES ON THE BODY 


In Figure 3, the input data listed an interpolation Mach number of 1.15, which is 
greater than the last body surface Mach number, 


Figure 7 shows the sonic line and stagnation streamline data derived by interpola- 
tion. Along the stagnation line, the local flow angle should be zero, except at the last 
point (the body stagnation point), where a singularity exists. The results shown in 


Figure 7 represent reasonable values. 


Figure 8 shows the resulting initial value line obtained by interpolation in the 
blunt body data. The initial message is caused by the fact that the interpolation Mach 
number was specified to be larger than the maximum Mach number on the body sur- 
face, calculated in the Blunt Body portion of the program. When this occurs, the 
program automatically resets the interpolation Mach number to be less than the final 


body surface Mach number, 


The listing of the freestream properties and geometric properties of the blunt 


body solution is followed by the interpolated drag coefficient, the table in which the 


interpolation was carried out, and the value used for interpolation. 


SONIC POINT OATA 


x R MACH NO VELOCITY TOT PRES PRESSURE TEMPERATURE Density OE, TA THETA 
© 120033€+00 = =.535332F +00 1,.0000010 44232226 +00 .796869F-02 .7661936902 .137400E202 .699786E*03 44.09680 67,20578 
eTS5O9IE+00 = .570K06F 600 2.0000009 .423722E +00 = .777957TF-02 .7HBOIOF 02 .137400E+02 ,699787E*03 45.80217 
e208209E +00 .604705F +00 1.0000038 .42)222E+00 .TS01TTF-02 .7I168SE*02 .137400E%02 .6997B87E+03 46.91614% 
2?6H752E 600 .633135F*00 1.000004% 421223600 .683925F-02 .657598E +02 .137400E*02 .699786E +03 47.51266 
2327921E*00 .6468910F 200 1.0000117 .62222HE 00 .613318F-02 .5A9TOSE*O2 .137400E*02 .699784F 203 67.97760 
2 3TSG4BE 600 .667356E200 1.0000374 .421736E100 .550363F-02 .S291T0ESO2 .137399ESN2 .699779E 403 48.24784 


STAGNATION POINT DATA 


x R MACH NO vELOCITy TOT PRES PRESSURE TEMPERATURE DENSITY OECTA THETA 
eLOSOHRE-06 .ETEO2Z2ZE-15 = 5. 3902817 = 1734664E 600 =.550365F-02 .9ORPOLECO2 .252942E +02 .102407E*0% -.06733 90,00000 
o262717E-01 0. 23068528 .137435F+00 .5$0365F-02 .938346E*02 .1553466E*02 .105354E*06 -.07108 
«SORRGOE-01 0. 2170013 .9RO25IE-O1 .550364F-02 .96931KHF 902 .158022E%02 .10TAZTE*O% -.14715 
e82GAP2F-0)  PPIOZ2F=125 — e BPASE2H «= 5607546E=-03 .550364F-02 .99IBZ0E*02 .161064E*02 .109K09E 206 9-2 44152 
2133946736200 ~.1D1022E-15 20013815 .A3TIAIE-03 = .550364F-02 -10016AF 203 .164834E 02 .110286E 204 42.6628 
eDI72L0F 600 ~.821626F-05 =. 00043RF =.4350G1E-03 .550363F-02 .1001696*03 .164AG9E*02 .110387E*04 BR.90918 


FIGURE 7. SONIC LINE AND STAGNATION STREAMLINE 


, e INTERPOLATION MACH NO, DECREASED 19 1.090 QUE TO LIMITATION ON BODY MACH NO. 


INITIAL VALUE LINE DATA 
MACH NO, GAMMA 
-8800000E°01 + 14000006 +01 


DELTA RBODY BODY 
«© 1372834€ +00 + 1000000E+01 + 1000000E +01 


coz 2622474E +00 BASED ON THE FOLLOWING TABLE 
yr) CPRIT) 
+682594E +00 »62569B8E +00 
+662436E +00 +601610E+00 
+641935E+00 +576398E +00 


4 INTERPOLATED FOR Y= 2579999E 600 


TVL DATAsSTARTING AT SHOCK 
x 


2 « Y DELTA MACH NO. PTs/eTO 
4 «2623376E +00 ~8140287E +00 240 74A896E +02 «1530358E +01 21297562F-01 
© 2835066E +00 «7T963341E +00 24143026E +02 +1421652E+01 2 1083439F-01 
. + 3082310E*00 »7713309F +00 242305 76E 602 o1318021E+01 290646206 -02 


 < 2 33TOSISE*O0 2 7445 109F +00 24348350E +02 ©1222897E +0) 761969) -02 
4 «3705602 +00 «7133409F +00 + GSIAOSGE 602 «114381 7E*0) 26422311 F-02 
+4051535E +00 +6 7T99995E +00 471 1418E 602 2 1089626E +01 2550 3634F-02 


a3 MASS-ENTROPY TABLE 

ae Oo. 255036 346-02 o1298246€-01 25504835E-02 2 3894721F -0) 255166536 ~02 
06491218E-01 ©5533738F-02 +9087705F-01 ©5562789E-02 2 1168619F +00 25601 756F-02 
+ 142R0468E +00 +5650840E-02 + 16P7717E +00 «S71L0297TE=02 «19473656 +00 «57804396 ~02 
«2207014E +00 2586)636F-02 224666636 +00 ©5954321E-02 «2726311F 200 «6088995E~02 
©2985 960E +00 o6176229F -02 «32456096 +00 +6306672E-02 «= 3505258F +00 2645105 7E-02 

‘ + 3764906E +00 26610207F-02 24026555E +00 267A5065E-02 +42R4206F +00 © + 6976605F ~02 
©4546 1B52E +00 + 7186039F -02 «480 3501F 200 © 7614636E=02 ~50631S0F +00 «766 3A26F ~02 
532279AF +00 © 7935213F-02 2552 2447F 600 28230579E=02 258420966 +00 oAS51913F-02 
+6101745E +00 +A901436F-02 263613936 +00 «92B1628E-02 246210426 00 29695 262E ~02 
e6BA0691F +00 © 10145446F-01 #7149339E +00 +1063566E=01 «73999RRE 200 ~11L16981F-01 
«76596 37F +00 o1175230F-01 « 79} 92A5F ONO «123880BE~0) «81460287F 600 21297562F-01 


FIGURE 8. INITIAL VALUE LINE 
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The final table is the mass-entropy table which is a listing of the shock point or- 
dinates (which define the swallowed mass) and the associated total pressures. Each 
line consists of three pairs of data, with each pair being a shock ordinate and total 
pressure. The first pair represents the stagnation streamline and the final pair the 


initial value line shock point. 


Figure 9 shows the first page of the Rotationally Symmetric Method of Character- 
istics output. The body geometric description occupies the upper portion of the out- 
put. The body is now located in a coordinate system centered at the body stagnation 


point. The blunt segment equation is then 


Y= 2Rx - BX 


where the values of R and B are given in the output. The second section is a special- 


ization of the general cubic equation 


Y = A(X-X,)° + B(X-X,)? + C(X-X,) +D 


BLUNT SEGMFNT 1 X2= ,826352F*00 Y2= .984R08F +00 xe= 41001646 +00 
R= .100000F*01 A= .100000€+01 

CURIC SEGMENT 2 x2= ,200000€+62 Y2= .2)9)75F +0) Kez .106326F 00 
SLOPE 1= 10.00000 DEG SLOPE 2= 10.00000 DEG 


= 0. C= 1763276 +00 D= .9R4B0BE +00 
POSAHOHCASEOHHHOREHOEEEDEDHEOEOHHYUOHOSOSOBONOTATIONALLY SYMMETRIC FLOWS SPCC OPEC ROOD HE EECO HEHE SHOOT HOODS ODOSOOROTOOOROOS 


33 ENTROPY VALUES READ 
==-=C=SUMMARY NO. Dee eee ee ee ee ee eee eee rrr 
x Y PT/PTO 
0233277 0713361 45.167 1.1438 0937454 2642231E-02 
o25R262 «706202 45,282 1.1389 0873937 eo SOSR44E-02 
2283815 2697910 45.740 1.0985 «848606 ©550363£-02 


ORAG COEF .= .643534E+00 (DELTA M)7(M TOTALD= 2427081F-01 


POINT x PTseTO 


2199768 2764517 1.00615 © 761069E-02 
0227292 2738517 2950355 0693776E-02 
0251657 0732413 #890620 0642690E-02 
2276805 «726505 2856779 oSIDHASE -02 
23039571 0717626 44.14) «822257 #550363E-02 


DRAG COEF .= .666944E+00 (DELTA M)Z(M TOTALD= o11)710F-01 


=--=-C-SUMMAPRY NO, 

FONT . PT/PTO 
2170948 2771331 42.3046 1.05331 2906462E-02 
2197288 2767619 42.4658 1.00756 ©8224613€-02 
2223760 2763010 42.521 2956394 27497766 =02 
#247648  T5B269 42.338 2699161 2693122E-02 
626954R 2753393 42.609 -667618 26473966 =02 
#294566 2766760 42.780 820562 25975 73E=02 
2325773 2738526 42.396 2 7K1742 25503636 -02 


DRAG COEF,= .691241E+00 (DELTA M)Z(M TOTAL)D= +78) AB0F-03 


FIGURE 9. ROTATIONALLY SYMMETRIC METHOD OF CHARACTERISTICS 
OUTPUT 
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where A, B, C, D are given in the output. For each segment, X2, Y2 are the coordi- 
nates of the final point. For segment 1, these coordinates represent the juncture be- 


tween the sphere and the cone of slope A-a@ (body slope minus angle of attack). 


Next is printed a reminder that this is rotationally symmetric flow, followed by the 
number of mass-entropy points read. This should correspond to the number prepared 


by the Initial Value Line portion of the program. 


The data for each characteristic line are grouped under C-SUMMARY headings 
which are numbered consecutively. The first few C-Summaries will start from points 


on the IVL and these will not have shock points. Eventually, however, the C-Summaries 


will have shock points as Point 1. This is denoted by the value of the shock angle 


being printed out as THETA SHOCK 1 = XXX just beneath the C-Summary number 
and by a 1 being printed at the extreme right of the shock point data. 


Each C-Summary contains, as a final item, the summed wave drag coefficient 


and the mass flow error. After the first few C-Summaries, this mass flow error 


should be much less than 0.10. Otherwise, a smaller step size is needed. 


ANGLE OF ATTACK =15,00 CONE ANGLE = 25.00000 DATA PLANE = 25773817 


6.800000 1.400000 15.000000 .57738174 
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FIGURE 10. 
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2.1932926 


MACH 
1.3527225 
1.4051956 
1,5603617 
1,7605251 
1.9743206 
2.1299215 
2.1854259 


MACH 
1,5055606 
1,.5014025 
1.7¢57297 
1.9429379 
2.16275646 
2.3€467610 
2.382A647 


MACH 
V.713R772 
1.7/69999 
1,9601956 
202 342H65, 
2.5009086 
2.693 7K77 
267620758 


MOCH 
1.996323¢4 
2.0764296 
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3.1780746 
3.26BR019A 


PT/PTO 
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#9083405 -02 
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+1256621€-01 
©1283465€-01 


PT/P TO 
©1275552€-01 
013462631€-01 
2155243356 -0) 
01926272E-01 
«23659766 -01 
o272A776E-01 
©2RAS18356-01 


PI/PTO 
27106778 -01 
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+? KIAIIBE-O1 
a 19A27 7G -Ot 
56100736 -01 
©70904196-01 
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THETA PSI 
26 o FPPGGP79=90.. OODDDD0D 
26,01400608-45, 34219559 
27,96334537 -5,194690644 
2890455539 29,03199324 
27.96334537 56.54005111 
26201400608 76.36295185 
24.99999979 96.00000000 


THETA PSI 
23.46334930-90.00000000 
25..065368 70-466 ,634661771 
28.26905526 -5,54567256 
3137634769 26,06426598 
32.75191677 51.81681664 
33.00740706 72.44252829 
32297651530 90.00000000 


THETA PSI 
22.28056799-90.00000000 
24.17716680-43,95265369 
28,09339206 -5, 34510087 
32207705509 25.310772%96 
34 662268828 50.32562417 
35. 7ASOTI3L 71.28336665 
36.10582939 90.00000000 


THETA Ps! 
21.57700346+90.00000000 
23674117593-43,5916588) 
28,36124020 -5,64989508 
32477665704 24.59163262 
36403404669 49,29875027 
37476813957 70.54981016 
38430402320 90.00000000 


THETA Pst 
21,11600824-90.00000000 
23.40411353-463, 30339389 
28459447790 -S.91061171 
33,58825526 23,79562061 
37,17560737 48,52)84246 
394331460276 70,017107R4 
40.00518452 90,00000000 


INITIAL VALUE SURFACE 


When the drag coefficient and the mass flow error no longer change, the C-Summaries 


no longer terminate on the body but out in the field. 


Figure 10 shows the final output, the initial value surface, derived by interpolation 
in the RSMoC-generated data. The data are presented "ring-wise," with each of the 
points being located on one of the specified meridian planes in the wind axis system. 
These data are also written on TAPE4, which can be defined as the system punch 
tape or a utility tape, for later reading by the the Three-Dimensional Method of 
Characteristics Program. 
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SECTION III 


OPERATIONAL ASPECTS : 


1, CORE AND TIME REQUIREMENTS 


The program requires approximately 100,000, words of memory to load and exe- 


cute on the CDC 6600 computer, 


The program requires 30-40 seconds of Central Processor Unit (CPU) time to 
compile. A normal run will require 15-20 seconds CPU time for execution. Very 
little difference has been found between the execution times generated by the OPT = 1 


and OPT = 2 compiler options. 


2. ERROR MESSAGES 


The built-in error messages in the program are listed Jelow. They are arranged 
i according to the OVERLAY segment and the program or subroutine in which they are 


generated. Probable causes and possible corrective measures are outlined. 


However, no such discussion can be complete. The program is too complex to 
allow every eventuality to be foreseen. When failures do occur in the Blunt Body 
portion of the program, it is probably due to a poor choice in the number of steps to 
the body. If the failure is in the Rotaticnally Symmetric Method of Characteristics 
portion of the program, a slight change in the step size, the number of points on the 
initial value line, or Mp ody for the 1VL will usually allow the program to continue. 


“4 In the latter case, a plot of the characterisitc network will usually reveal the problem. 


A failure in the Rotationally Symmetric Method of Characteristics portion of the 
program will result in the print option, INCASE, being reset to a value of 4 and the 
: control being set back to repeat the calculation of the present characteristic. The 
physical appearance is the sudden onset of extremely detailed print, If the failure 
occurs in the region governed by the initial value line (since the program cannot re- 
wind the input tape) the next point on the IVL is read in, and a detailed print of that 
calculation up to but not including the body point is given. By plotting the characteristic 


data, the troublesome point is easily spotted. 


| 
f 


OVERLAY (BLUNT, 1, 0) 
Program BLUNTS 


COMPUTED FEWER THAN 3 BODY POINTS TWICE 


Not enough shock points. Increase number of shock points. 
Too many steps to the body. Decrease number of steps. 


FAILED TO FIND A SUPERSONIC BODY POINT IN 4 TRIES 
Too few shock points. Increase number of shock points. 
FAILED IN G CONVERGENCE 
No solution for the n-step multiplier. Change number of steps to the body 


EXCEEDED 9 ITERATIONS FOR BODY 


Convergence criterion too small. Relax convergence criterion. 
Roughness in the relationship between Bhody 294 Bshock: Plot the individual 
data, and establish a starting value for Bghock.- 


ZERO OR NEGATIVE G 


Too many steps to the body. Decrease number of steps. 


Subroutine BODY 
NOG ROOT FOR DELTA 


Poorly defined body due to roughness and/or few body points. Increase number 
of points on the shock and/or decrease number of steps to the body. 
INSUFFICIENT NUMBER OF BODY POINTS FOR CURVE FIT 


’ Not enough body points. Increase number of shock points and/or decrease num- 
ber of steps to the body. 


Subroutine IDTHER 


4 ROOT2 FAILURE, TR= _BM= AM= 


Calculated static pressure is greater than calculated total pressure. Arrange 
number of points on the shock and number of steps to the body so that AN is 
approximately equal to A€ 


e 


OVERLAY (BLUNT, 2,0) 
Program IVLS 
OID NOT FIND SUPERSONIC POINT M= __ 


Not enough supersonic flow region in the Blunt Body portion of the program. 
Increase number of points along the shock. 
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Subroutine BBOUT 


EXTRAPOLATED SONIC POINT 


Not enough supersonic flow region in the Blunt Body portion of the program. 
Decrease number of steps to the body. 


OVERLAY (BLUNT, 3,0) 
Subroutine BODI 


X1 LIES BEYOND LAST SEGMENT 
Probable error in body description and IVL geometry. Check input data 


XB= 


POINT LIES BEYOND END OF BODY X3= 
Body description geometry is not sufficient for given data points. Plot character~ 
istics to check for geometric incompatibilities, and check input data. 

ITERATION FAILED 


No convergence on body point. Plot each iteration and check for errors. 


IM EXCEEDED 25 
Body point Mach number is less than 1.0. Increase Mpody in input data. 


Subroutine FIELD 


EXCEEDED 3*#50 ITERATIONS IN FYELD POINT 
Probable geometric inconsistency. Check input data, plot characteristics data. 


FAILED TO CONVERGE S3 IN SPECIAL LOOP 
Poor definition of the relationship between the stream function and total pressure. 
Check mass-entropy table, may be a read error by computer. 

FAILED TO CONVERGE D3 IN SLOW LOOP 


bis’ Geometric inconsistency. Check input data. 
a 
; DR FAILURE IN FL(ELD 
First approximation gave the correct answer, making the residual almost zero. 
The subsequent iteration was therefore thrown off. Change step size slightly. 


FAILED IN FIELD» IMM= __ P3= 


Field point Mach number is less than 1.0. Plot characteristics data, check for 
geometric inconsistencies. ; 


TRIED 10 TIMES TO CONVERGE POINT BUT STILL SUBSONIC 
Field point Mach number is less than 1.0. Increase Mbody in input data. 
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**ENOTE**** CONVERGENCE CRITERIA MULTIPLIED BY 10** , FOR POINT 


Poor convergence at field point. Plot characteristics data to determine 
troublesome point. 


OSCILLATION IN FIELD FOR DELS 


Convergence problem in determining the local value of the stream function. 
Examine mass-entropy table for smoothness; may be a read error by computer. 


Subroutine ENTROP 


AFTER LAST ENTRY IN S/R TABLE. GMASS= 2SR SET TO 


LAST ENTRY,» »W/CORRESPONDING MASS 


»SR SET TO 


BEFORE FIRST ENTRY IN S/R TABLE. GMASS= 
FIRST ENTRY» »W/CORRESPONDING MASS —____ 


« GMASS SET TO 


BEFORE FIRST ENTRY IN S/R TABLE. S/R= 
FIRST ENTRY,» eW/CORRES. S/R= 


- GMASS SET TO 


BEYOND LAST ENTRY IN S/R TABLE. S/R= 
LAST ENTRY» »W/CORRES. S/R= 


These statements are advisory in nature, indicating that the mass~entropy table 
and the local calculations are poorly matched. While the calculation will! continue, 
the mass-entropy table, the geometric parameters, and the input data should be 
examined for errors. 


Subroutine WAVE 


FAILED TO CONVERGE SHOCK ANGLE IN 50 TRIES 


The subroutine was given an incorrect deflection angle across the shock. Trace 
origin of the deflection angle, using check-out print. 


4 Subroutine SHOCK 


SM IS ZERO OR NEGATIVE IN RSHOCK 


A negative step size is specified or complications arising from the mesh build-up 
resulted in a negative step size. Check input data and/or plot characteristics 
data. 


OSCILLATED THROUGH NGATE THEN IGATE 


OSCILLATED THROUGH IGATE THEN NGATE 


Both of the above messages indicate that base point data for the shock point 

solution are being determined by extrapolation rather than interpolation. Plot 
the progress of the shock point solution to find out where the incorrect extra- 
polation occurs. 


gr Ty i. se TEM ak ARTI 


EXCEEDED 50 ITERATIONS FOR POINT 4 


Geometric inconsistencies in the mesh build-up. Plot characteristic mesh to 
locate trouble spot. ig 


EXCEEDED 75 ITERATIONS FOR SHOCK 


Poor initial data. Plot characteristic mesh to locate problem area. 
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SECTION IV 


LOGICAL STRUCTURE 


1. OVERLAY STRUCTURE, READ/WRITE LOGIC 


The logic of the program is controlled by the master program and uses four first 
level OVERLAY segments. The master program consists of a main program to 
call each OVERLAY segment and one general use subroutine. The Blunt Body Program, 
the Initial Value Line Interpolation Program, the Rotationally Symmetric Method of 
Characteristics Program, and the Initial Value Surface Interpolation Program 
each comprise one OVERLAY segment. Figure 11 shows the OVERLAY structure 


and the storage requirements. 


Figure 12 shows the read/write logic structure. TAPE2 and TAPES would be the 
same device if minor changes were made in the program logic. TAPE4 is shown as 
a dotted output device to indicate that this could be either the system punch file or 


a scratch tape, as the user specifies. 


2. INTERDEPENDENCE OF SUBROUTINES 


The Calling-Called matrices for three of the four OVERLAY segments are shown 
in Figure 13, The segment comprised of the Initial Value Surface Program consists 
of one main program which in turn calls two general purpose subroutines. This 
matrix, due to its simplicity, is not shown in Figure 13. 

3. LISTING 


The complete listing of the IVS program is given below, following Figure 13, 


LY 


U] 
ee oe ee . . og eee. 
LOM ce meta ENCANA J, Ci OR ASN ae ie ms 


Overlay (BLUNT, 0, 0) 


Program [VALUE 


Overlay (BLUNT. 1. 0) | Overlay (BLUNT, 2. 0)} Overlay (BLUNT, 3, 0){ Overlay (BLUNT, 4, 0) 


Program Program 
BLUNTS IVLS 


FIGURE 11. OVERLAY STRUCTURE 


FIGURE 12. READ/WRITE LOGIC 
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hee ee 


a 


(c) Program MOCS 


FIGURE 13. CALLING-CALLED MATRIX FOR IVS PROGRAM 
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SECTION V 


FUNCTIONAL DESCRIPTION OF SUBROUTINES 


OVERLAY (BLUNT, 0,0) 


IVALUE - Main program. Calls each subordinate program in sequence. 


ENTRPG - General Lagrangian interpolation. 


OVERLAY (BLUNT, 1, 0) 


BLUNTS - Main program for blunt body solution. Reads in data, sets up initial 
parameters, and governs course of iteration for the correct shock 
wave shape. 


BODY - Determines equation for body, given the set of body points, and de- 
termines the coefficients for the equations describing the Mach 
number variation along the body. 


XY - Calculates the Cartesian coordinates of a field, body or shock point. 


DIFATE - Performs a seven point central difference calculation to determine 
the first and second derivatives of data along a given shock, field 
or body line. 


IDTHER - Determines the ideal gas thermodynamic properties of a given point 
in the flow field. 


’ IMAIN - Carries out the inverse blunt body solution based on the parameters 
specified by BLUNTS 


OVERLAY (BLUNT, 2, 0) 


4 IVLS - Main program for determining the properties along a line through 
: the blunt body flow field in order to start the rotationally symmetric 
‘ method of characteristics solution 
DETER - Interpolates for all properties at a given point. 
i BBOUT - Prints out the final correct blunt body solution. 


OVERLAY (BLUNT, 3,0) 


7 MOCS - Main program for rotationally symmetric method of characteristics 
solution. 
BODI - Carries out body point solution. 


| 


OPER8 - Organizes the field and body point solutions and determines when the 
calculations have been carried sufficiently far downstream. 

FIELD - Carries out field point solution. 

VMACH - Calculates Mach number from static and total pressures. 

ENTROP - Keeps track of the relationship between total pressure and stream 
function, and, given one of these, determines the other. 

TSON - Calculates the shock wave angle necessary to develop a downstream 
Mach number of unity. 

DIF - Performs linear interpolation along characteristics. 

CHANGE ~ Promotes and demotes characteristics during the calculation. 

INPUT ~ Reads the input data and makes the necessary dimensional changes. 

SHAPE ~ Generates the values of the coefficients defining the body and the 
step sizes required. 

DROP - Checks the geometric relationship between points on the initial value 
line and just downstream of the bow shock and drops any unnecessary 
points. 

OUTPUT - Arranges and prints the program output. 

WAVE - Solves the Rankine-Hugoniot relations to determine the properties 
downstream of a shock wave. 

: SHOCK - Carries out shock point solution. 

MESH - Controls step size along the bow shock in order to achieve desired 

step size along the body. 

PAGE - Keeps track of lines printed and skips to a new page when necessary. 
Also prints page number in upper right hand corner. 

‘ ANGLE - Computes Mach angle from Mach number. 

VELOC - Calculates V/Viyax from Mach number. 


OVERLAY (BLUNT, 4, 0) 


‘4 IVSES - Main program for the interpolation of the initial value surface. 
; ROTATE - Rotates a vector through a specified angle. 

by 

; SECOND - A specialized second order Lagrangian interpolation. 


124 


( 


PART 2: THE THREE-DIMENSIONAL METHOD OF 
CHARACTERISTICS PROGRAM 


SECTION I 


INTRODUCTION 


The Three-Dimensional Method of Characteristics Program calculates inviscid 
supersonic flows about smooth three-dimensional bodies at angles of attack. Given 
an initial value surface (IVS), freestream conditions, and a properly described 
body, the program determines the flow field downstream of the IVS until a user- 
specified station is reached. The program can also run back-to-back with the 
Initial Value Surface Program. In this case, the IVS is automatically provided by the 
Initial Value Surface Program, and only the freestream conditions and the body de- 


scription need be provided by the user. 


The program uses two temporary files TAPE2 and TAPE4, and on options it reads 
the IVS from TAPE4 and writes the output on TAPE1, 


{ The Three-Dimensional Method of Character:*#ics Program has computed a vari- 
ety of flow fields and has proven to be efficient and versatile. Certain limitations of 
the program are noted here. The local Mach number must be greater than 1,0 every- 
where. Hence the freestream must be supersonic, and the configuration must be 


such that no subsonic region exists at the juncture between the canopy and the fuselage 


c or between the wing leading edge and fuselage. The body must be smooth, without 
- any surface slope discontinuity (i.e., at any point on the body a unique normal exists). 
4 Different numbers of data points per ring (up to a maximum of 48 points) may be 


assigned in different body regions (up to 20 regions). 


rr 


SECTION II 


BODY DESCRIPTION 


Every configuration has a number of generating lines, such as the upper profile, 
the lower profile, the maximum breadth line, or the wing leading edge. In the pre- 
sent body description procedure, each generating line is divided into a number of 
segments to permit each segment to be described by a conic-section curve. At 
each cross section of the configuration, simple analytic curves, e.g., the ellipse or 
cubic, connect any two adjacent generating lines to form the contour of the corss 
section. The configuration is thus described analytically by simple low-order curves. 
For a smooth body a unique normal to the surface exists everywhere, and this con- 
dition usually requires slope continuity at the junctures between two contour curves 
or two segments of a generating line. 

The fuselage is located in a right-handed coordinate system where the Y-axis 
is aligned with the fuselage axis; the X-axis is spanwise and the Z-axis is up. All 
generating lines are represented by a general curve fit of conic sections in several 
segments. The conic-section curve takes the form 


Zz 9 1/2 
(3) = PY+Q+s¢(RY + SY + r) 


where SG=+1, A straight line is a special case with R=S=T=0. Each curve can be 
divided into as many segments as necessary to provide adequate body description, Each 
segment must be continuous with the previous segment and with very few exceptions the 
slope must be continuous at the junctures to satisfy the requirement of a unique normal 
to the surface, 


A typical cross section of the wing-body configuration is shown in Figure 14, 
Through each point marked by a dot, a generating line which is a function of Y 
passes. For input cards, the following Table shows the curve number correspond- 


ing to the two projections of each of the generating lines. 
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Segment Type 


Co Ellipse 

Z om " 
GI " 
nl " 
L¥F¢ = mt 
wae. A 

; S04. Circle 

BG: ut 
ie Cubic 


Straight Line 
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Designation Designation 


For instance, curve number 3 denotes the projection on the Y~Z plane of the generating 
line that passes through point L for both Phase I and Phase II designations. Aside 
from the wing, the only difference between these two phases is in curve 7, Whilea 
shape factor b was used in Phase I to alter the lower fuselage, a side flat section 


~ NNNNN 
Qa | a Oo Fig 


NN MM 


ow 


was used in Phase II. 


As shown in Figure 14, the contour of the cross section begins with a straight line 
representing the canopy flat. The canopy contour from points C to ¢ is circular but 
can be elliptic in general. The upper fuselage is represented by an elliptic curve 
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from points U to M, A straight side flat from points M to Fe joins the upper fuselage 


f toL, A 


straight line from point L to the centerline describes the bottom flat. For the 


to the lower fuselage, which is also represented by an elliptic curve from F 


description of the wing, straight lines EG and FH represent theupper and lower sur- 
faces of the wing, respectively. Partial ellipses GI and HI complete the wing des- 
cription near the leading edge I. The wing-body blending is effected by circular 
ares 34 and 56 with radii 


ry and "> respectively, 


In general, the configuration are divided into a number of sections since some 
of the generating lines can only be described by a composite of conic-section curves. 
Some of the conic-section curves, however, may remain valid throughout a few 
sections. To avoid repetition of the input data, the body description subroutine 
FIGURE reads the data in the following manner (see also Section IV, Input/Output 
Description). For the first body section, the fist card specifies the end of the 
section and seven conic-section curves. This is followed by 7 cards of body descrip- 
tion coefficients, each describing one curve. For a subsequent body section, the 
first card specifies the end of that particular section and the total number, KT, of 
conic-section curves that require new values of the coefficients for that section. 


Only KT cards are then used to supply these new values. 


In Figure 14, the projections of lines passing through E, F, G, H and I on the 
Z-Y plane and X-Y plane are input quantities for defining the wing. In order for the 
partial ellipses GI and HI to exist, points G and H must be located within a certain 
range, which depends on the relative position of these five points. When the wing span 
is very small, it is difficult to input both projections of lines through G and H such 
that these points are located within acceptable ranges. In such cases, the X coordi- 
nates of G and H are calculated inside the program to satisfy the range requirement. 
As shown in the sketch, 


Z 


E 
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when the projections Xp Ze X Zr and Zo are given through input, XG can be 
calculated when the ratio R = Gb/ab is specified. 


Indeed, 
RGb ab... “G2 
=f0° ee, eee 
2.2. (x. — 
be Pees en oe tp 
and X, = X;-Gb = X,+ —S ZT. 


when fitting the generating lines, a fuselage station is selected ahead of which the 
projection Xo is not fitted because XGq will be generated by the program according 
to the above formular with a constant R taken at this fuselage station and input to 
the program. Since at this station the derivative of XG with respect to Y can be 
calculated from the formula, this derivative should be used to fit the projection Xo 
for fuselage stations immediately behind the selected station. 
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SECTION III 


GRID ARRANGEMENT 


Since the geometry of the wing-body configuration is fairly complicated, a pro- 
cedure was developed to arrange the grid points for efficient calculation. The con- 
tour of a wing-body cross section has a number of regions that undergo sharp 
changes of curvature, where a very fine grid structure is needed. On the other hand, 


in regions only mild changes of curvature occur, a coarse grid structure suffices. 


A sketch of a typical contour is shown in Figure 15, together with the grid points 
and the labeling of controlling key points. Quantities K)) HS ireiclie’s K are specified 


through input data; for instance, if K, = 20, it means K, is the 20th point counting 


from the lower profile. Quantities Ky to Kio are sbieiees by adding to K, the 
differences between K, and Ke to K., respectively; for instance, if Ke = 17, then 

Ky = 23. The quantity Kg specifies the number of points to be added in regions near 
the wing-body juncture and the wing leading edge; for instance, Kg = 2 as illustrated. 
541? Key to 
Ky, Kyo-1 to Kj 5 not counting the additional points specified by Ky. Anywhere 
else double spacings are assigned. The grid size below the wing leading edge, where 


Basically, single grid spacings are assigned from K, to K,, Ky to K 


K, is located, is in general different from that above the leading edge and the difference 
depends on the value of K,. Additional points specified by Kg are added to each single 
grid space from K, to K,, K, to Ky, Kio to Kyi: Care and some experience are needed 
in assigning the K's, However, the program has the capability to center the series of 
points from K 4 t0 K; (and hence Kyo to K, » around the wing-body juncture. The 

group of points from K, to K, is designed to cover the fuselage-canopy juncture. All 
K's are set to zero at fuselage stations before the canopy. Between fuselage-canopy 


junction and the wing-body juncture only K, and K, are assigned proper values while 


all other K's are set to zero, 


a’ 


e 
@ 


FIGURE 15, GRID STRUCTURE AND CONTROLLING DATA POINTS 
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SECTION IV 


INPUT/OUTPUT DESCRIPTION 


1, INPUT DESCRIPTION 


The Three-Dimensional Method of Characteristics Program offers some input and 
output options. The program can run using the IVS data on punched cards or it can run 
back-to-back with the initial Value Surface Program, which provides the IVS data 
through a temporary file, TAPE4. 


The input data cards required by the Three-Dimensional Method of Characteristics 


Program are shown in Figure 16, They are described in detail below. 


Card 
No Variable Format Description 
@ E LIT 18A4 Problem identification 
aS 
DP F10.5 Convergence criterion for P, 0,1 to 10 
DTHETA F10.5 Convergence criterion for &* 0.1 to 10 
DPSI F10.5 Convergence criterion for ¥* 0,1 to 10° 
EPSLN F10.5 Convergence criterion for A, 0.1 to 107° 
@ DYC F10.5 Initial step size which satisfies the Courant- 
Friedrichs-Lewy stability condition, i.e., 
DYC-tanys<J2 d, where As is the local Mach 
angle and d is the shortest distance between 
two data points. If DYC=0, the program will 
generate an initial step size. 
YSTOP F10.5 Computation stops when Y 2 YSTOP. YSTOP 
cannot be set equal to 0. 
STEPM F10.5 Maximum step size (when STEM > 0) 
wincstt F10,5 Fuselage station at which the wing begins (110 


for the sample case) 


*The relationship between the flow angle © or W and the coordinate system is shown in 
Figure 16a, 


tVariables or cards so marked are not needed for the Phase I program. 
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Card 
No Variable Format Description 


NSURF 15 Total number of surfaces desired, counting 
from IVS, even for a restart from a data 
surface 


LINE 15 Number of data rings per surface. LINE is 
updated by data from card 7 for a restart from 
cards 


NPC 15 Number of data points per ring. NPC is updated 
by data from card 7 for a restart from cards 


NR 15 Number of regions having different NPC. 
1s NR $ 20 


IFT 15 Control for reading IVS data. 


= 0, IVS data is on cards 
If1FT 4° 1s data is on TAPES 


Also, if IFT=K< 0, then number of rings per 
data surface is increased by {K| during restart 
from TAPE1. |K|s 2 is recommended. 


ISTART 15 Restart control* 


<0, restart from cards 
; If ISTART , =0, initial run 
| =1I>0, restart from Ith data 
©} surface on TAPE] 


KP4 15 Output control. If KP4=K > 0, every data sur- 
face is written on TAPE], with every Kth sur- 
face printed out. If KP4=K #0, every data 
surface is printed out. If K=0, every 10th 
data surface is written on TAPE]. If K< 0, 
every Kth data surface is punched out on cards, 


¢ NS 15 Identifying number of the surface on which 
print level KPRINT (1) will apply 


ie NL 15 Identifying number of the ring on data surface 
4 NS along which print level KPRINT (2) will 
apply 
NP 6) Identifying number of the data point along ring 


NL on data surface NS at which print level 
KPRINT (3) will apply 


KPRINT (1) 15 Print level control usually 0 


( =0, normal print 
If KPRINT (1) , =1, moderate default print 
= 2, extensive default print 


KPRINT (2) 15 Print level control, usually 0 or 1 
KPRINT (3) 15 Print level control. See Section V 


*For restart, insure ISTART #0, or data on TAPE] will be ruined 
j **See Section V for warning on print level 
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i per ring to be calculated is specified by the 
BPSET(Gl-9) a variable NPSET(9I-8), and the grid size is 
NPSET (91-7) 15 specified by NPSET(9I-7) . . . NPSET(91) 

: aa Section II] where Ky denotes NPSET (91-7) 


Ky denotes NPSET (90))* 


Card 

No Variable Format Description 

YSET(D F10.5 For all Y s YSET (1), the number of points 
epeerent 15 -. 


@ There now follow cards, I = 2, NR 


FMSTR F10.5 Freestream Mach number 
GF (1) F10.5 Ratio of specific heats 
© ALFA F10.5 Angle of attack, in degrees 
YO F10.5 Y-Coordinate of initial value surface. 
Needed only for an initial run using IVS 
data from cards. May be left blank for 
all other cases, 
' YA F10.5 Y-coordinate of center of spherical nose 
ZA F10.5 Z-coordinate of center of spherical nose 
SCALE F10.5 Scale factor (radius of spherical nose) 
© WINGA F10.5 Wing area, Sq. Ft. (reference area for 
; Cpw) 
DBLNL F10.5 Distance between points 5 and F in Figure 14 
(0.5 for sample case) 
¢ 
5. + NCQ 15 Number of sections needed to fit the 
us @ blending curve (#22 for the sample case) 
4 by conic sections 
J 12 Section number 
+ CQ(1, J) E15. 7 Fuselage station representing lower limit 
® of the section 
CQ(2, J) E15.7 Fuselage station representing higher limit 


of the section 


*Note that for the Phase I program only Ky and K, are used. 


tVariables or cards so marked are not seoted for the Phase I program, 
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Card Variable Format Description 
No 
IC 12 Curve number for blending (#22 for the 
sample case) 
@' CQ(3, J) E15,7 Conic coefficient P(J) 
CQ(7, J) E15.7 Conic coefficient T(J) 
CQ(8, J) F2 Conic coefficient SG(J) 
| 
® There now follow sets of cards, J = 2, NCQ, each set describing a section. | 
@ If ISTART 20, omit cards 10, 11, 12, 13 
@® Note: Cards 10, 11, 12, 13 are required only for a restart using data on 


punched cards. They are generated by the 3DMoC Program when the flow 
field data is punched out through the KP4 option. Value of the restart 
parameters are also part of the printout for each data surface. The variable 
INCRES may be set at a different value as explained below. 


ISURF 15 
LINE 15 
i ILL 15 Restart parameters 
: KURE 15 
MESS 15 
INCRES 15 If LINE is the number of rings per data 
surface during a restart from cards, then 
‘ the new number will be LINE + INCRES, 
INCRES € 2 is recommended, 
NPC 15 
. KI(1) 15 Restart parameters* 
nf K1(s)t 5 
‘f SINE E15. 8 
j COSINE E15.8 Restart parameters 
® $2 E15.8 
C2 E15, 8 


*When making a restart after an initial ran, be sure to remove IVS data cards, if 
they were present for the initial run. 
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Card Variable Format Description 

No 
DYC E15. 8 
STEP 

@) re BtGve Restart Parameters 
YO E15, 8 
DRAG E15. 8 

©) ALIFT E15.8 
TORQUE E15. 8 Restart Parameters 


@ GO TO card 15 unless EE = o-one SET =), \or 


——=————_— 


ISTART <0 
A(J,3,) E15.8 Z-coordinate of data point 
A(J,1,D E15.8 X-coordinate of data point 
©) A(J,7,D F10.7 Mach number at data point 
A(J,8, J) E15.8 Total pressure ratio, Pye 
THETT(J,D F12.8 Flow angle 8 
PSIT (J, I) F13.8 Flow angle ¥ 


@ There now follow cards, J = 2, NPC which describe a ring on the data 
surface 


@ There now follow sets of cards, I = 2, LINE with each set describing a 
ring on the data surface 


YMAX F10.5 Y-coordinate of the end of the first body 
4 section described by the coefficients to 
; follow. 
© N 15 Number of body description curves for 
bg (©) the body sections. Must be set equal to 7. 
: mop? 15 = 1, fuselage only; = 2, fuselage - canopy; 
hy’ = 8, fuselage - wing; = 4 fuselage - canopy- 
wing 
RATIO? F10.5 Ratio for partial ellipse (see Section I) 


If Ratio > 0, curves 19 and 20 are generated 
by the program. ** 


**For the sample case, RATIO > 0 between F.S. 110 and F.S. 220 
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Card 
No 


; iD) 


Variable 


L 


P(L) 
Q(L) 
R(L) 
S(L) 
T(L) 
SG(L) 


@ There now follow cards, L = 2, 6 


Format 


I2 


E15.8 
E15.8 
E15. 8 
E15. 8 
E15. 8 
F2.0 


Description 


Identifying number of the body description 
curve 


Body description coefficients (See Section II) 


@ If the body is described by only one section, no further data is required 


YMAX 


P(M) 
Q(M) 
R(M) 
S(M) 
T(M) 
SG(M) 


F10.5 


15 


15 
F10.5 


12 


E15.8 
E15.8 
E15. 8 
E15. 8 
E15, 8 
F2.0 


Y-coordinate of the end of the next body 
section described by the coefficients to 
follow 


Total number of body description curves 
that need new values of P, Q, R, S, T, SG 
for this body section 


See card 15 
See card 15 


Identifying number of the body description 
curve which requires new coefficients 


Body description coefficients 


@ There now follow a total of KT-1 cards, giving the new coefficients for this 


body section 


@ There now follow sets of cards, starting with card type 17, each set 
describing the subsequent body sections to the end of the configuration 


Ea apt. ce Sl ck RIT Re 


2. OUTPUT DESCRIPTION 


A number of output options are available through the control index KP4, which is 
part of the input data. If KP4- K >0 , every data surface is written on TAPE] with 
every Kth data surface printed out for sampling. If KP4 =K <0, every data surface is 
printed out, while retaining certain data surfaces for a restart (see also Section V, 
Operational Aspects). If K = 0, every 10th data surface is written on TAPE]. If K is 


negative, then every Kth surface is punched on cards. 


When any output is written on TAPE], the user should either use a physical tape 
as TAPE], through a request control card, or copy TAPE] to a physical tape 
as many times as there are data surfaces through a COPYCF control card before 
termination. Otherwise, information on TAPE], which is a scratch disc file, will be 
lost upon termination of the run. A restart at a later time can be made from the physical 
tape. The data surfaces on the tape can be printed out by copying them to OUTPUT 
through a control card COPYCF (NAME, OUTPUT, N), where NAME is the tape designa- 
tion and N is the number of data surfaces to be printed out. The tape can also be taken 


to a peripheral processor for printing of the data surfaces. 


Figure 17 shows a normal printout for a data surface. X, Y, Z, are the coordinates 
of a data point (see Figure 15 for the coordinate system). P is the static pressure and 
PT is the total pressure, with the subscript 0 referring to freestream quantities. The 
angles OUTWASH and UPWASH, in degrees, are included since they are more meaning- 
ful than & and W. These angles are related to the direction cosines u, v, w of the 
velocity vector by: OUTWASH = atan(u/v), UPWASH = atan(w/v), @= acos(v), and 
Y > atan(w/u). 


Figure 18 shows a data surface where the CFL stability condition was violated at 
certain data points. When the velocity vector is aligned with the Y-axis, certain flow 
angles are undefined. When this condition is detected by the program, the coordinate 
system is automatically rotated and the particular data surface is recomputed, An ex- 


ample of such occurrences is shown in Figure 19. 


When the normal output is written on TAPE1, every Kth surface is printed out for 
sampling. An example of this is shown inFigure 20, where the 5th data surface is 
printed out completely, while the upper and lower centerline pressure coefficients on 


the body surface are printed out for data surfaces 6 through 9. 
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DATA SURFACE NO 50 


*e*TO USE THIS AS AN INITIAL VALUE SURFACE» KEY PUNCH FIRST 6 COLUMNS OF DATA AS INPUT CARDS 


SURFACE NO SOe 


es eee 


~e263273191 F601 
= 22493021 7E eG? 
=. 209186936+01 
-el47622286 +0} 
«7266751 0F +00 
2 l0587877E +00 
+ F10R9550E +00 
+1628) 831E +0) 
e21667R7T2E°01 
225) 36693E +6) 
e2O323191F +01 


es eens 
22797T6743E +01 
= 22652468126 +01 
=e 223462966 01 
= 015948933€ +01 
2799410416200 
275389643E-01 
210540670F +01 
2193969536 +01 
2266571226 +01 
23101281 7E +01 
2 3258 746B9E +01 


ey ee ees 
=~ e294665306E +01 
~028122352E 01 
~o2369T403E 201 
-016971259E +01 
2961701266 *00 
©63255148E-01 
© 12081033€+01 
022666129E +01 
eo 3E260502F +01 
© 368891 06E +01 
@ SARGAGODE OR 


ee ee 
BIISI4SIE*G? 
29713223E+01 
250238 70E +0} 
el TIIIITSF+Ol 
—e PL IG7HBLE $00 
2641561036 -01 
o137)364)6+0) 
2 26015709E +01 
#361590 74E +0) 
2422829 36F $01 
2452394236 201 


ecJorccose 


33755736E +01 
= 232009045E +01 
270031790! 
19631LA37E*OL 
=0999STIIOF *00 
oS53P17HBF-01 
«© 15389353E+0) 
© 29624946F +01 
e169 131F 691 
o49621380F 401 
oS22117H0E°01 


TO USE AS IVSe+ 
50 s 6 

¢. 

219505687€*00 


STATION Y = 


Br en Ore 
=.20385704E-11 
eBoGhI7RIE O00 
.15977245t +01 
221793025E°0) 
2532A4 1E 0) 
226 30077IE +01 
2 24696906E +01 
«20T6I46) FESO) 
+ 14945006E+01 
+ 7TRI2Z31465E +00 
= .95486356E=12 


aan Rectan tara 
~.JO7TOROTTE=11 
2903464 39E+00 
17137388E +01 
2 235474670E+01 
o2772T064E +01 
o2938727TIESO1 
2 2A469322E+01 
226631095E+01 
oT TASTIGSE+O1 
©9424 3RA0E +00 
-, 1A669713E~11 


SsusweKece=s= 
- 30161 T6SE-11 
297196 IRRE 200 
eo LAS6SOTIIECO! 
2 254444620E +01 
3020011 7E+01 
«32468202E +01 
2 32270727E +01 
o2AOTSIOIE SOL 
020 764533E 20} 
o110405467E+01 
=, 228877 3eE et) 


Semoun Xacosee 
=. 30659667E-)) 
2 106250GRE +01 
2) 9BO995)E +0) 
o276099H0E +01 
© 3276942ARE +0) 
2 356090466E +01 
2 359552616 +01 
o J17846251E+01 
223736922E +01 
©12700193E+01 
= e2390A9K3E-11 


KeSesehvcssas 
= .3033553RE-11 
011265813E+01 
221450391E+01 
2298062466E +01 
« ISAGAROIE*O1 
239350575E +01 
+40362166E +01 
2 36054105E +01 
227127502E +01 
o146580109E +01 
~e2272967T0E-11 


MaChH NO, 

1,6777867 
1,7229238 
1,6123940 
129760670 
2218861846 
22430667) 
2.6926111 
2.9180800 
320737217 
3215274146 
321740096 


MACH NO, 

26656985 
226963417 
2.7828730 
2.9377430 
3,1232295 
323255777 
35724237 
37451265 
38498956 
328991642 
349115300 


MACH NO. 

362299571 
3e2B10111 
Bot 31 PROT 
326658618 
30 9GREL 94 
42601023 
4 oSS61987 
4.7705603 
4 9044262 
49667372 
4, 9836407 


MACH NO. 

324733950 
325470307 
347633794 
421003025 
424957342 
4,8916615 
5.2997R05 
5.9966917 
5.7875031 
5.591096 
529236606 


MACH NO. 

326197886 
367110667 
369770099 
438746689 
428511827 
522953026 
567528929 
620922807 
663029156 
6.9300267 
664719569 


269602655E +01 


BOOY POINTS 


---PT/PI0---- 
26703995 -02 
«47039950 -02 
67039950F -02 
247039950F -02 
+67039950F -02 
26793995S0F -02 
267039950F -02 
267039950F -02 
267039950F -02 
267039950F -02 
267039950F -02 


THE TAs DFG 

14,.99999998 
35, 360214658 
16,19329874 
17.2 7TLLT287 
17,880337§81 
18.121359S6 
17.93319395 
17,187874680 
16,17656692 
15.32772 32) 
14,99999998 


FIELD POINTS 


==-PT/PT0~--- 
27292557E-01 
227043901F-01 
026363194F-01 
«25123766 -01 
22333604RE-01 
©21341877F-01 
«196336266 -61 
© 17649843E-01 
«16066152 -01 
«15053060F-01 
2 146693739F-O01 


THE TAs DEG 

13,31310132 
13.5189278)1 
14,04846027 
14, 71087442 
15, 34037100 
15.78929729 
16.048065649 
16.065346175 
1584616070 
15,58465241 
15,4754044¢ 


FIELD POINTS 


---PT/PTO~--=- 
oS4EREGSOF-Ot 
2SST4RGTOE-01 
25653697 IF =O} 
258 36AR923F-01 
@59265227E-02 
2585240336-01 
«SSSITI3I“F-0) 
+5099ROS6F -01 
24732) 3275-0) 
246204069F-01 
«431 34643F-01 


THETAs DFG 
B1.92 969348 
12.13770859 
12.40522A87)} 
13.2710718) 
13.A40631)6 
14.31354213 
16.919692346 
15.45497728 
3574299546 
15.84236371 
15, 85967158 


FIELD POINTS 


---PT/PT0---~ 
+68) 79800F-01 
+ 7059A330F-01 
»77823781F-01 
2 88867509F -01 
elO1L1L6011F +00 
211230892F +00 
«11960076¢ +00 
«12265760€ +00 
«1229462AF +00 
«12282711 +00 
©12123839F +00 


THETA, DFG 

10.47699753 
10.78636306 
11.51956? 16 
12.35586470 
13.17638674 
13,.87056906 
14,90520463 
15,93036388 
16.72504136 
17.23176789 
17,40628270 


SHOCK POINTS 


---PT/PTO---~ THETAs DEG 
2714676644E-01 8,39261285 
eTSTS3O0TF-01 8.90781S22 


2892032366 -01 
211296150F 00 
21462361AF +00 
2178721295 +00 
«218824676 +00 
«251470 79E 00 
2272921 36€ +00 
«28622905 +00 
+29076667F +00 


10.10634147 
11.49100404 
12.90737422 
14.13077635 
15.86885527 
17,49489797 
18,71557603 
19,477287632 
19.73673985 


PUNCH THE 3 CARDS+ PLACE BEFORE ROUY PT CARDS 


Qo 0 6 0 @ n 
«10000000E+0!1 0. + 10000000E +01 
+50000000E-01 .69602455F +01 


FIGURE 17. NORMAL DATA SURFACE 


--PS1 s0EG-- 
-90.00000000 
-58,90000772 
29, 9661664) 
=4.23829123 
18.05198180 
3734813305 
54235713544 
6729081624) 
77.93828499 
8489024934 
90200000000 


=-PS1+DEG-- 
=90-00000000 
-61.23774330 
-34,16076513 
-9, 76328100 
12,6912862496 
30.206754653 
47664526664 
61.98GIL1L26 
73.40760919 
R2.40581379 
90,00000000 


~-PS1+DEG-- 
~90.00000000 
~60. 35033663 
33207388640 
Fe 8257228 
232333)1646 
28200018866 
45234500773 
59 .46656577 
TLeLIONT79L 
81.03921897 
90.00000000 


--PS1+DEG-~ 
-90.00000000 
-57.80054023 
29411237350 
-S.02107111 
14, 383462325 
29.43056499 
45.10801612 
57,97844 30% 
69.18231679 
79.71300680 
90.00000000 


--PS1+0EG-~ 
-90.00000000 
-53.85709166 
=23.804044642 
~ 46116161 
17,34666876 
30648329776 
44,64861526 
57,29123861 
68,47196915 
74,28847527 
90.00000000 


ween C Penn n = mee -P/P0- n= == =P /PT0---OUTWASHIIP WASH 
0.00-15.00 
8,06-13,22 
2 10SS4F-02 14,12 -8,75 
+834616F-03 17613 -1241 


- 33202€+00 
o31S®1E*00 
~27245E*00 
o21421E 900 
+ 157G4E*00 
»llO77E*00 
+ 76390E-01 
+55742E-01 
244925E-01 
240276t-01 
«S9LIOE-01 


oe --C0----+ 
© 339K0E #00 
2322946900 
22?B%2E +69 
2218126 *00 
«15945E*90 
11L277E°00 
«76228E-0) 
25S1S6E-01 
243773E-9) 
2 38294E-01 
2 I6690E-0) 


wece~Cpeer- 
23305E *00 
2314273 900 
o27T096E*00 
#21374E°00 
+1S7A0E*00 
211376E*00 
«77408E-01 
2SO7TS4E=9) 
245400E-01 
2 396°0F-01 
2 37B899E-01 


ee 
231228E *60 
+ 29604E +00 
«2562300 
«20325600 
215656E +00 
«11709€ 00 
1856936-41 
.67S2 76-01 
eSTIVIE-01 
SLB TSE-O1 
+501 706-01 


te 
«27652E*00 
22664 7E 400 
223297E*00 
+19310E*60 
o 1STRIE*OO 
et32P4E S00 
2109236*00 
2958S 7E-0) 
oABS4RE-01 
2 843°6E=-9) 
@BI00RE-O01 


2 62426E 602 
259428 +02 
eD1SO0TE +02 
#40632E +02 
+ 300546E +02 
2 21493E +02 
151336 +02 
11 II3E+02 
eFSLITESOL 
oR®4S16E +0) 
+82356E +01 


----P/P0---- 
243829E +02 
250 76BE +62 
2 52492E 6 02 
2413556 +02 
« 30500F +02 
221864E +02 
«TSVIGE 602 
o11206E +02 
«909R5E +0) 
2A0B48E +01 
«77RBIE +O) 


=---P/P0---- 
e623 37E +02 
2591 3HE 02 
eSLI31E +02 
240452F +02 
2 30195E +02 
222046E +02 
o1S321E +02 
#1 1S00E +02 
2 93995t 601 
eAIISTE +O 
SOLN7E +O} 


2---P/P0s--- 
2 SB7T6E 602 
0SS938F +02 
SRG 2al +G2 
« REZ IE +02 
2295935 07 
ao P2664E our 
« 16AS44 +02 
2 1 36956 G2 
oe LLA03E +02 
of OS97E +02 
« LO2R2E +07 


a3 -P PQ ene 
S2160E +02 
4 O9IOE 202 
2441026 02 
« JATZ6E 662 
+ 30200E+02 
«252826 +02 
2212286 +02 
« DBI35E 602 
1 7368E +02 
«16601F +02 
e 1GASTE Oe 


21 2AL6E-02 
122006 -02 


617006 -03 
#441256 -03 
+ 3J1068F -63 
+23225F-03 
eG 7E-03 
o17351F-03 
© 16908F -03 


17.05 
14.58 
10.68 
6.64 
3447 
1.40 
6-00 


5.71 
1).23 
14.76 
15.99 
15.A6 
15.27 
15.00 


=P/PT0-~-OUT WASH IP WASH 
0200-1343) 
6.60-11.90 


et3t0“t~02 
«124726 ~02 
«10776F ~02 
~84901F-03 
«626)61 -03 
+44BATE~03 
+33 029F ~03 
+23003£-03 
»18679F-03 
+ 165985 -03 
» LSIROE -03 


1).70 
14.51 
15,04 
13.73 
10.97 
7.70 
4.63 
2.11 
0.00 


76.00 
72.55 
3.18 
6.10 
12.00 
14.26 
15.22 
15.45 
1S .4R 


=P/PT0---OUT WASHIID WASH 


ol2757t -02 
121461F-02 
+L 0497E-02 
© 83048F -03 
+61991F-03 
+45261F 03 
2 31454F-03 
236106 -03 
2 19297F-03 
eI TIL 36-03 
« 16468F -03 


0.00- 


11.92 


6007-16659 


10.66 
1361) 
13.59 
12.79 
10.61 
£.00 
Sel 
2253 
6.00 


ADD 
2016 
2.7? 
6.63 
10.73 
13.49 
14.99 
15.66 
15,46 


USPTO --- OUT wA SHI WAS: 


ol 2067t -02 
oe L16R4E-02 
2996) 3-03 
« TIARO -G3 
200753 -03 
465284 -63 
« SGAG1F $03 
2277056 <0) 
«298217 203 
2b ISHF =09 
«211 09F -03 


=P/PT0---0UT was 


= 10706F -02 
«10251 -02 
«905416 -02 
eTSIVTEH03 
+62000F-03 
253903 -02 
2S 3ISROF-03 
2 38425 =03 
«35656 -03 
2 S40R2F -03 
© 3 3SR2F-03 


0.00- 


5.Re 
10.19 
beet 
12.78 
te.t« 
10.54% 
Aa) 
©.10 
3037 
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SECTION V 


OPERATIONAL ASPECTS 


1. CORE AND TIME REQUIREMENTS 


The program requires a core size of 60,000, to compile, 1000,000, to load, and 
61, 300, to execute, all on the CDC 6600 computer. 


The program requires approximately 30 seconds of Central Processor Unit (CPU) 
time to compile on the FTN compiler, using OPT = 1. A typical run will require 
approximately 10 minutes of CPU time for the fuselage of Phase I and one hour for the 
wing-body of Phase II. There is practically no difference in the execution times 
generated by the OPT = 1 and OPT = 2 compiler options. Therefore, OPT = 1 is 


recommended because of its somewhat faster compilation. 


2. START OPTIONS 


a. Normal Start 


Provisions have been made to run both the Initial Value Surface Program and 
the Three-Dimensional Method of Characteristics Program in one pass. This is done 
by combining the control cards; putting the Initial Value Surface Program and data 
before the Three-Dimensional Method of Characteristics Program and data, and setting 
ISTART = 0 and IFT =1. 


During execution the 3DMoC program will automatically read the initial value 
; surface from disc file TAPE4 where it was placed by the IVS program. It is also 
possible to run the IVS program and the 3DMoC program separately, using a physical 
tape to transfer the IVS information. In this case the IVS data on TAPE4 are recorded 
on a physical tape through a COPYCF card in the Initial Value Surface Program 
at the end of the computation. In the Three-Dimensional Method of Characteristics 
Program, set IFT = 1 and ISTART = 0 and transfer the IVS data from the physical 
tape to disc file TAPE4 through a COPYCF card before computation begins. The 
physical tape (TAPE1) on which all data surfaces will be written, can double up 

| 


for the job of transferring the IVS data. 
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b. Restart 


Normal output can either be written on a physical tape (TAPE1) or be printed out. 
When a restart is made from the Nth data surface on the tape, set ISTART =N. Be 
sure that ISTART #0, otherwise the data stored on the tape will be ruined. 


When the normal output is printed out, a restart can be made by one of the following 
methods. When the output control KP4 is set to zero, every 10th data surface will be 
written on TAPE}, which can then be used for a restart by setting ISTART = m(10). 

A restart can be made from any completed data surface by key-punching the data 

cards from the listing and making a restart from cards with ISTART set equal to 

any negative integer. In key-punching the data, ignore any line with literal information; 
simply key~punch the first six columns of data and stack the last four cards first, as 
described in Section II]. During restart, the number of data rings per data surface 

can be increased. Set INCRES = K for restarting from cards or set IFT = -K for 
restarting from tape. The number of data rings N will be increased to N+K. It is 


recommended to limit K to 2 or less. 


3. DEFAULT PRINT 
A printout option is built into the program to print intermediate results at a 
particular point, along a particular data ring or on a particular data surface at the 
lower level KPRINT = 1 or at the higher level KPRINT =2. This may be done for 
check-out purposes when other means have failed to reveal a source of computational 
} trouble. Level 2 gives extensive default print, while level 1 gives a moderate amount 
of default print. As a rule, up to level 2 is used at a data point, while level 1 is 


c used along a data ring. Otherwise, the amount of printout obtained could be extensive. 


In case a failure occurs during a restart from a physical tape, it is useful to keep 
track of the data that are stored on the tape. For this purpose, KPRINT(3) can be 
set equal to 3 to print out the first 40 columns of the stored data as they are read in 
during the restart. Other intermediate results will be printed out at the higher level 2 


at point NP, along ring NL, on data surface NS when KPRINT(3) is set to 3. 


4. ERROR MESSAGES 


If the program runs all the way, it is a fairly good assumption that the results are 
good. Usually, if something goes wrong, the computation will stop sooner or later 
with an error message. Default print can be obtained, if necessary, to assist in 
locating the source of the error. In the case that the normal output is on TAPE], 
and the computation stops at a data surface which is not printed out, it will be necessary 
to restart from the previous data surface and set KP4 = 1 to have the partial data sur- 
face printed out. This will locate the troublesome data point and detailed default print 


may then be obtained for that point. 


Error messages, unusual stops and other program messages are summarized below 
below. They are grouped according to the program or subroutine in which they are 


generated, 


Program FUSLAG 


ILLNATURED PSI AND DELTA CASE NO—, ROTATE COORDINATES AND 
RECOMPUTE SURFACE 


Whenever the local velocity is aligned with the Y-axis, the angles W and 6 are 
undefined. When this happens, the coordinate system is rotated around the 
Z-axis through an angle of 0.1 radian. The above message is printed out (for 
the Ith such occurrence) and the Nth data surface is recomputed. 


DATA PCINTS WILL BE RE-SPACED FOR SURFACE | 


This message will appear when (1) the distribution of data points along a body 
ring becomes noticeably uneven, (2) a data ring is automatically added by the 
program, (3) a data point computation violates the CFL stability condition 


WRONG IVS NO —_, NO —_ SPECIFIED OR FAIL END FILE | 


When a restart is made from TAPE1, this message indicates either the program 
is unable to locate and read from the specified data surface or the arrangement of 
data on TAPE] is unusual so that an END OF FILE was encountered. If this 
happens, the contents of TAPE] should be printed out. Alternatively, when 
KPRINT (3) is set to 3, the first 40 columns of each line information on TAPE] 
will be printed out as the program FUSLAG is reading TAPE]. 


Subroutine BULK 


BULK POINT HAS FAILED TO CONVERGE AFTER 25 ITERATIONS 


Computation is stopped. 
This happens very rarely. A decrease in the magnitude of the convergence 
criteria may remedy the problem. 
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Subroutine BASEPT 


BASE PCINT __ VIOLATES CFL 
This advisory message appears whenever the CFL stability condition is 
violated at basepoint I. The number following the message is the factor 
L (see the discussion on basepoint location in Volume f. 


INTERSECTION FAILURE, BASEPOINT NO __ 


Computation is stopped. 

This message indicates that the Mach conoid fails to properly intersect the 

line joining basepoint I to the center data point. The first three numbers 

following the message are coefficients A, B, C, of the quadratic equation 

(see discussion of basepoint location in Volume J) and the last number is 
B?-AC. When this happens at the juncture point of two body sections 

cheek the continuity of the body description. If the body description is new, 

check FIGURE, paying attention to the formula and coding of the body normal. 

In high gradient areas, this message usually means more data points are 

required, 


MACH CONE CONVERGENCE FAILURE IN BASEPOINT NO __ 
Computation is stopped. 
This seldom happens. A relaxing of the convergence criterion on Mach angle 
(increase in the value of EPSLN) helps. 

BASE PT __ EXTRAPOLATES», STEP DECREASED 


The number following this message is the factor L (see discussion on basepoint 
location in Volume I). If this message appears in a series prior to an abnormal 
termination, an increase in the number of data rings (and in some cases, in the 
number of data points) may help. 


Subroutine COM PAT 


CONVERGENCE FAILURE IN COMPAT 


Computation is stopped. 

This has never happened. If it should occur, try to increase the number of 
iterations by modifying the source deck. If it should persist, obtain default 
print for analysis of the problem. 


Subroutine GMT RY 


CONVERGENCE FAELURE IN GMTRY 


Computation is stopped. 
Check the body description and subroutine FIGURE, 
Subroutine FIGURE 


$$$$% YOU HAVE GONE ALL THE WAY--FEELING NICE, I BET $$$$$ 


This message indieates that the goal has been accomplished 
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SECTION VI 


LOGICAL STRUCTURE 


1, INTERDEPENDENCE OF SUBROUTINES 


The Calling-Called matrix for the program is shown in Figure 21. 


2. LISTINGS 


Two complete listings of the Three-Dimensional Method of Characteristics 
Program are given below, following Figure 21. The first listing is for the computer 
program developed in Phase I for the calculation of flow fields over a fuselage and 
canopy configuration. The second listing is for the computer program developed in 
Phase II for the calculation of flow fields about a fuselage-canopy-wing configuration. 
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SECTION VII 


FUNCTIONAL DESCRIPTION OF SUBROUTINES 


Coordinates the three basic modes of computation: body point, field 

point, and shock point calculations. It establishes the position of a 

new point at each iteration and checks to see if the convergence cri- 
teria are satisfied. 


Picks the appropriate data points from the last computed data surface 
for each of the three basic modes of computation. 


Determines the location of the base point (the intersection of a line 
joining two data points with the backward facing Mach conoid) and the 
flow properties at this point. 


Computes the values of the coefficients of the compatibility relation 
in the generalized finite-difference form for each of the three basic 
computations. 


Solves the nonlinear equation resulting from manipulation of the com- 
patibility relations for the shock point computation. 


Finds a new body point and the unit normal at that point. 


Describes a slab-delta-wing type body through six conic-section 
curves and provides the unit normal at any specified body point. 


Transforms to and from a body-normal oriented coordinate system. 
Rotates a two-dimensional vector through an angle. 
Solves three simultaneous linear equations. 


Calculates the average between the locally computed shock normal 
and the shock normal determined from the shape of the shock surface 
and recomputes the local velocity behind the shock, using the average 
normal, 


Sets the mirror point properties across the plane of symmetry. 
Normalizes a vector. 

Performs all output options. 

Rearranges the data points along the data rings. 

Rearranges the data rings on a data surface. 

Provides the coefficients of a cubic spline for data interpolation. 
Performs data interpolation through a cubie spline. 

Rotates the X-Y corrdinates. 


Computes Y and Y' from conic-section coefficients. 


Computes th: cubic contour near the fuselage-canopy juncture. 
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ELIPSE 
WINGF 
PCRCL 
PELIPS 
CONPIC 
XSUBG 
AERC 
SETCP 
SECOND 


Computes the elliptic contour of the fuselage. 

Computes the upper or lower wing contour. 

Computes the upper blending circular contour, 

Computes the partial ellipse contour near the wing leading edge. 
Computes the upper limit point of the upper blending contour. 
Determines the X abscissa of point G of Figure 15. 

Computes the pressure coefficient and moments. 

Collects relevent data for AERO. 


Makes interpolation using the 2nd order Lagrange's formula. 
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PART 3: THE SURFACE FIT AND MINIMUM SEARCH PROGRAM 


SECTION I 


INTRODUCTION 


Latin Squares are used to sample configurations for wave drag calculations. 
Corresponding to each cell of the Square, a configuration exists which can be 
represented by a set of body describing input data, Both 5 x 5 and 3 x 3 Latin 
Squares are used: the former is for Phase I study and the latter, Phase II study. 

, can be treated by a 5 x 5 Latin Square and four variables 


6 
Lyreecer dy by a3 x 3 Latin Square. These reduced variables z, are related to 


Six variables Zyreeeee sd 


the physical variables Xx through the formula 


xX; See Z. (Xs -Xj_. 

_ ~ Imax Imin i| Imax i itin) 

X. = ee sf 
1 2 n-2 


where n = 6 for a5 x 5 Latin Square and 4 for a 3 x 3 Latin Square. The advantage 
of the reduced variables is that a program written in terms of them is valid for any 


physical variables with any ranges of variation. 


When the wave drag coefficients Cow have been calculated as input data, 
the Surface Fit and Minimum Search program is used to obtain the coefficients 


by bis and bi of the wave drag equation* 


Cc = bo + >», z+ Dy 252, 


Dw 


through a least wquare procedure; in other words, the program is used to fit a 


hypersurface in an n dimensional space, When the type of constraints is specified, 
4 the program is used to yield a set of values for the Zs that corresponds to a minimum 


wave drag body subject to the given constraints. 


‘ *i and j range from 1 ton: however, as explained in Volume I, not every possible 
| term is included. 


a 


SECTION II 


GUIDE TO THE SEARCH PROCEDURE 


In the search procedure, wave drag coefficients are computed using the wave 
drag equation at points covering the entire region bounded by the given ranges of 
variables in the n-dimensional space. The point which has the lease wave drag is taken 
to define the minimum wave drag configuration, The accuracy of locating the mini- 
mum is therefore dependent on the resolution of the numerical network. Accurate 
location of the minimum can be obtained efficiently by a search-by-steps technique 
where the search is first conducted using coarser grids to determine an approximate 
location, and then repeated in a smaller region centered around the location of the 
previously obtained minimum. This process can be repeated as many times as 
needed to obtain the desired accuracy. The major advantage is that the given con- 


straints can be satisfied easily by rejecting points which violate the constraints. 
1. SEARCH STEPS 


The searching procedure is summarized as follows: 


a. Select the ranges of search by specifying the lower and upper limits of 


each of the n number of reduced variables:* 


{ n = 4 for 3 x 3 Latin Square 
n = 6 for 5x 5 Latin Square 


b. Select the resolution for each of the variables by specifying the numbers 
of grid points N, (j = 1, n) of that variable. The resolution is given by 


Z\ -(z 
ie u), ( L), 
resolution - —— AY 
J 
* (21) = O(J), (24) = E(J), and 25 = NP( J) in the input data described in 
Section III. 
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The total number of points for each search cycle is equal to the product 


T 
of Nj} i es 


total points = II N; 
=i 


For example, if 10 points are specified for all variables, the total number 


of points for n = 6 is 10°, To keep the searching time within practical 
limit, Ri must therefore be selected niacin 4 so that the total number of 


ee points is within reasonable limits * Using the data of a. and b., 


the program computes the n_ tables of variables: 


( i=, N. 


lia), - 


(74), 


so that all configurations defined by every possible combination of the 
n variables, 


can be checked for the given constraints (e.g., the volume). The con- 
figuration which has the minimum wave drag coefficient computed by the 


wave drag equation is considered the minimum wave drag configuration. 


If the search is to be repeated, select the new ranges of the n variables 


by redefining the lower and upper limits as, 


Ve = zZ!-A,, (z 
1), J J u), 


where a are the variables for the minimum wave drag configuration deter - 
mined previously, and A, are arbitrary increments of Z selected to define a 
smaller region of search, In situations where (2,) or (Zu) defined by 

i , 


*See Section IV for time estimates. 


above relations are beyond the ranges of consideration, they are then set 


equal to the minimum or maximum of the range. 


d. Repeat the procedure from step b, until the desired accuracy is obtained. 


2. TYPES OF CONSTRAINTS 


Four types of geometric constraints are included in order to satisfy most 
design needs. These are: 


a. One or more of the six geometric variables take assigned values within 
the ranges of variation. For instance, the length can be kept the same 
as that of the baseline or the width can take 90% of the baseline width. 


b. At one or more fuselage stations the configuration cross sections contain 
geometric curves prescribed by tabulated data of width X versus elevation Z. 
This type of constraint is useful for the placement of equipment. In the 
case of the radome constraint, only the radome radius and the elevation 


of the center need to be specified. 


c. At one or more fuselage stations the configuration cross sections satisfy 
minimum area requirements. The minimum area can be either a given 


area or a certain percentage of the baseline cross-sectional area. 


d. The configuration satisfies a minimum volume requirement between two 
given stations. A given volume or a certain percentage of the baseline 


volume can be imposed. 
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SECTION III 


INPUT/OUTPUT DESCRIPTION 


Two different sizes of Latin Squares are used in Phase I and Phase II; 
the input and output of these two phases are described in this section. 


The input data cards for Phase I program LATIN1 and Phase II program LATIN2 


| 1. INPUT DESCRIPTION 
are shown in Figures 22 and 23. The cards are described below: 
/ 
, 
| 
| 


Phase I Program LATIN1 Input 


Card Variable Format Description 

No 

gE 

> @ SYM(1) 3 x A2 | Alphanumeric symbols to identify reduced 

. : ( variables Zyr ss ees Ze in print out 

SYM(6) 3x 2A 

Q) IOR 15 Number of terms in the Wave Drag equation 

7 IPRINT 15 Print out control 

| 

j 

be ( = 0, Minimum Print 

IF IPRINT ) =1, detail print about 

| Be surface fit 

ay 

: (3) KIND 15 Constraint Control 

, 

, - = 0, no constraint; or con- 
straints imposed on one 
or more of the variables 

= 1, local area-ratio constraint 
If KIND 


= 2, local tabulated X-vrs-Z 
| ; constraint 


= 3, randome constraint 


= 4, volume constraint 


@ Use card 4, if KIND = 1 or 2 
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Card 


© 


Variable 


FS 


RATIO 


Format 


F9. 2 


F952 


NEALE T CTRL EN TE LIS BE PEE BA fot pe ss enn SENS 


Desc ription 


Fuselage station at which cross-sectional 
area constraint is imposed 


Ratio of desired area to baseline cross- 
sectional area at this F.S, if0< RATIO s 2 


Desired cross-sectional area at this F.S. 
if RATIO > 2 


Not used if KIND = 2; may be left blank | 


® Add cards 4, and 59 if KIND = 2 


NPFS 


Z.FS(1) 
XFS(1) 


ZFS(NPFS 
XFS(NP FS) 


15 


B9C2 
¥F9. 2 
F9, 2 
F9.2 


@Use card 4. if KIND = 3 


FS 


FC 
RD 


@Use cards 4, and 5 


4 
FS] 


FS2 


NFSI 


FS(2) 


FS(NFSI-1) 


Blank Card 


F9. 2 


F9, 2 
F9, 2 


FLO. 2 


F10. 2 


F10, 2 


15 


| 


—_— 


Number of points in the table 


Elevation Z(1) of the first point 
Width X(1) of the first point 


Elevation Z(NPFS) of the NPFS*® point 
Width X(NPFS) of the NPFst® point 


Fuselage station at which the radome con- 
straint is imposed 


Z value of radome center 


Radius of Radome 


4 if KIND = 4 


Fuselage stations between which the volume 
constraints is imposed. Note that if FS1 < 50 
the complete nose is included. 


Ratio of desired volume to baseline volume 
if 0< ANY ¥s2 
Desired volume if ANY >2 


Number of specified fuselage stations at 
which cross-sectional areas are calculated 
for volume integration 


Fuselage Stations at which cross-sectional 
areas are calculated for volume integration 


Must be present to terminate the input 
string 


Card Variable Format Description 


No 
@~ J 12 Cell number of the Latin Square 
ANY F10.0 Wave Drag Coefficient for Jth cell 
e@There follow 24 cards like card 7 
i I2 Identification number (starting from 26) 
for additional input of wave drag coefficient 
ANY F10.0 Wave drag coefficient for Jth addition 
(8) V(1) F10.0 Value of reduced variable z, for Jth 
: : addition 
V(6) F10.0 Value of reduced variable Ze for Jth 
addition 
® There follow additional cards like card 8 
(3) Blank card Must be present to terminate the input 
string. 
elf KIND 4 4, go to card 13 
J 12 Cell number of the Latin Square 
ANY F10.0 Nose volume up to F.S. 50 for Jth cell* 
e@ There follow 24 cards like card 10 
J 12 Identification number for additional nose 
@) volume input 
ANY F10.0 Nose volume for Jth addition 
@ There follow additional cards like card 11 
@2) Blank card Must be present to terminate the input 
string 
NP(1) 15 Number of evenly spaced grid points for 
@) $ : reduced variable Zy for the search process 
NP(6) 15 Number of evenly spaced grid points for 
reduced variable Ze for the search process 
@) O(1) F9,. 2 Lower limit of the search interval for zy 
\4 ° ° 
O16) FO, 2, Lower limit of the search interval for Ze 
Seetow IV for volume calculations, 
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Card Variable Format Desc ription 
No 
© E(1) F9.2 Upper limit of the search interval for Zy 
E(6) 9. 2 Upper limit of the search interval for Ze 


Phase II Program LATIN2 Input 


Card Variable Format Description 
No 

SYM(1 

@) : (1) $ ie ( Alphanumeric symbols to identify reduced 
syM(4) Fe he ( variables Zyy +e y Zin print out 
IOR 15 Number of terms in the wave drag 

(2) equation 
IPRINT 15 Print out control 


= 0, minimum print 


If IPRINT = 1, detail print about 
surface fit 


KINDT 15 Constraint control 


/ 


= 0, no constraints; or con- 


straints imposed on one 
@) If KINDT | or more of the variables 


= 1, volume constraint 


VOLBIN F10.0 Volume increment desired (i.e., the 
desired volume minus the baseline volume)* 


> 


@ J I2 Cell number of the Latin Square 
NY F10.0 Wave drag coefficient for Jth cell 


e There follow 8 cards like card 4 


J I2 Identification number (starting from 10) 
G) for additional input of wave drag coefficients 
ANY F10.0 Wave drag coefficient for Jth addition 


*See Section IV for volume calculations. 
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Variable Format Description 
V(1) F10.0 Value of reduced variable zy for Jth 
: ° addition 
V (4) F10.0 Value of reduced variable Z4 for Jth 
addition 


e There follow additional cards like card 5 


Blank Card Must be present to terminate the input 
string. 


elf KINDT =0, go to card 10 


J 12 Cell number of the Latin Square 


ANY F10.0 Volume increment for Jth cell (i.e., 
volume of Jth configuration minus the 
baseline volume) 


eThere follow 8 cards like card 7 


J 12 Identification number for additional 
volume increment input 
ANY F10.0 Volume increment for Jth addition 


eThere follow additional cards like card 8 


Blank card Must be present to terminate the input 
string 
NP(1) 15 Number of evenly spaced grid points for 
2 Fé reduced variable Z4 for the search process 


NP (4) 15 Number of evenly spaced grid points for 
reduced variable Z4 for the search process 


AARON Anas rT A IN ow 


O(1) F9.2 Lower limit of the search interval for zy 
04) F9. 2 Lower limit of the search interval for 4 
E(1) F9,2 Upper limit of the search interval for zy 
E(4) F9.2 Upper limit of the search interval for z 


4 


ltl enitnth thscecttehe cates anions tnt iain inh asad ce ee owe . caro og 


hn 


OUTPUT DESCRIPTION 
Sample printouts of LATINI (5 x 5 Latin Square) and LATIN2 (3 x 3 Latin 
Square) are shown in Figures 24 through 27. Figure 24 shows the results of the 
drag coefficient surface fit for the 5 x 5 Latin Square. The input data shown in 
Figure 24 includes the drag coefficient for the basic 25 Latin Square cells plus 
7 additional configurations including the baseline (No. 26). With the 7 additional 
drag coefficients, 6 additional terms were added to the wave drag equation (see 
Appendix A of Volume 1), Figure 25 shows the 5 x 5 Latin Square surface fit of 
the nose volume. The printout is similar to that of the drag coefficient. Figure | 
26 shows the results of a minimum wave drag configuration search with volume 
constraint. Here the search was made keeping the variable z, at a constant value of 


1 


0, while the remaining 5 variables z - » 2, were varied through 5 discreet 


points between -2.0 and 2.0. That ‘a the resolution of each of the 5 searching 
variables in this search cycle is equal to 1.0. The minimum wave drag configura- 
tion, within the given ranges of variables and resolution is printed along with its 
cross-sectional areas and volume. The number of configurations which do not 
satisfy the volume constraint is indicated by the number of points rejected. Figure 
27 shows a sample output of the 3 x 3 Latin Square surface fit and the minimum 
value search results with volume constraint. The printout is similar to that of the 
5x 5 Latin Square results. Here, however, the fitted volume is the complete 

fuselage canopy and wing, and the numerical values shown are the differences from 
the baseline. 


(18) 
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(21) 
(22) 
#23) 


(ie 
( 6) 
(Iiy 
(16) 
(21) 


(1) 
( 6) 
all) 
(16) 
(21) 
(26) 
(31) 


[i 3) 
( 6) 
(ily 
(16) 
> (21) 
(26) 
#31) 


FIGURE 24, SURFACE FIT OF WAVE DRAG COEFFICIENTS (5 x 5 LATIN SQUARE) 


A 
VARTABLES ARE cistk6 
ORDER OF FITe « « +23 
CONSTRAINT = 4=VOLUME 


INPUT DATAscee 


ADDITIONAL TERMS ARE, 


reoperer 
POrirpr 


8.722590 
5.970170 
5.724920 
5.8d4690 
5952560 
7.625800 
7.254040 
7.105100 
5.439380 
>. 760620 
6.526840 
12. 6.673050 
13. 66591540 
16 6.861430 
1s 7.610590 
16 8.217730 
17 6.688300 
18 6.438170 
19 6.618100 
20 $.939600 
21 7.630460 


—ASORNFN FUN 


22 8.222750 


23 6.625220 
26 7.248460 
25 5.876170 
26 96 694370 
27 6.572600 
28° 6.095690 
29 22955290 
30° S.863430 
31 66814490 
32 6.118840 


*WAVE DRAGE 


AR 


“ “UM VALU EaRCH PROGRAM 
£ IN| Iv y te es ve 96 


COEFFICIENTS OF WAVE DRAG EQUATION 


6.589280 (2) =.3965690 
=.2050504E-01 ¢( 7) 2688001 
e2906077E-01 (he) =, 2480445E=-9) 
=.2082480E=-01 (17) -,4708812E-02 
~286279S5E-01 (22) ~1614850E-01 
STO O€v= .7IBLIE~01 
FITTED wAve DRAG 
8.657905 ¢ 2y 5.936778 
7.716832 € ib) 7.130032 
6.609727 (le) 6.476029 
8.263872 (17) 8.683026 
7.645751 (22) 8.136629 
64782620 (27) 6.607058 
6.801438 (32) 6.143416 
ERROR OF FIT 
~.646RG54E“01 ( 2) =-,3339190E-01 
«9103227E-01 ( 7y  =#,1246009%% 
«8288685E-01 (lé)  =.1964207 
246146162E-01 (17) =, 4473909E-02 
eIS29119E-01 (22) -,859210CE=01 
»8824958E-01 (27) «3445 787E=01 
-.1305239E=-01 (32) »24S57550E-01 


(3) 
( 8) 
413) 
(18) 
(23) 


(3) 
( 8) 
(13) 
(16) 
(23) 
(28) 


1952752 

2 3790653E-°01 
-5690S92E=-01 
4236248601 
© 1834393E-02 


S.756353 
7199410 
6.589280 
6.355805 
8.601087 
6.056381 


©3143257E-01 
FG IOFFIEMOL 
= e22599A6E-02 
= 8236504E-01 
= e2613316E°01 
-23930931E-0) 
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+ 206461 7E-02 
= 64525385E-03 
«1303170€-01 
1710669E-01 


5.968195 
5.445890 
6.947679 
6.736629 
7.299810 
5.948000 


©8350464E-01 
26510390€=02 
»8624890E-01 
21187293 
05134994E-01 
-27290010E-02 


«S) 
(10) 
as) 
(20) 


OS) 
(10) 
ais) 
(20) 
(25) 
(30) 


($) 
10) 
(1S) 
(20) 
(25? 
(30) 


4552590 

41 79657E-01 
-.146S3514E-01 
= .4901663€-01 


5.927520 
52702036 
72551826 
52841928 
5.939397 
58628546 


-.2503977E-01 
~.58586)0F-01 
-,5876602E-01 
= .9767153E-01 

«6322689E-0) 
-.5763780E-03 


ie 


INPUT DATA.,ee *VOLUME® 

27969.00 

24190.00 

22492.,00 

27633.00 

31198,00 

25953.00 

26918.00 

31499,00 

25273.00 

30202.00 

22460.00 

26212.00 

29430.00 

34228.00 

37811.00 

28259.00 

33819.00 

29035.00 

32614.00 

30919.,00 

20603.00 

31917,00 

35621.00 

36351.00 

31103.00 

25076.00 

25055.00 

23565.00 

21192.00 

20452.00 

27049.00 

: 271946,.00 

bp 29697209 ( 2) 1468.d06 (3) 854.0985 ( 4) +39.72363 ( S) 1778,2346 
( 6) 169,246} (7) 14232439 ( 8) -195.7796 «9 56.59233 (10) 202.7722 
(11) =184,8959 (le) 111,8339 (13) -189.0630 (14) 260.7751 (15) 108.9358 
(16) 230.2987 (7) 325,6316 (18) -546,28646 (19) 114.5688 (20) 288.7734 
(21) 372.5243 (22) -251 2878 (23) 320.3092 ( 
STO DEV= 212.15 
€ By 28199.48 23933.67 ( 3) 22326.73 (4) 278631.99 (5S) 31019.,06 
( 6) 25911469 ( 6) 31465.42 (9) 25708.52 (10) 30171.99 
ql) 22483.77 (13) 29697.09 (14) 33702.02 «isp 37972,46 
(16) 282346.93 33408.26 (18) 29127216 (19) 32664.66 (20) 30918.35 
(21) 26743674 31698.c4% (23) 35829.78 (26) 36640.07 (25) 30867,88 
(26) 25439263 25164.94 (26) 23504.00 (29) 21128.77 (30) 204446,52 
(31) 26999285 27234,34 


ERROR oF FIT 

(dl 230.4782 (2) 7256,3349 -165.2720 (4) 196.9875 -178.9608 
(6) -61,31367 ( 7) =56.65509 -33,57555 (9) 435.5181 =30,01195 
(11) 23.77196 (ley) =254,3200 267.0949 (14) -525.9835 161.6449 
(16) -26,07026 (17) =610,7430 92,16195 (19) $0,65513 = 6518263 
(21) 140, 7384 (22) -218,7605 208.7771 (26) 2890659 -235.1225 
(26)  363,6349 (27) 109, 9404 -61,00450 (29) -63,22808 7.476213 
(31) -69,14077 (32) 40.34173 


FIGURE 25, SURFACE FIT OF NOSE VOLUME (5 x 5 LATIN SQUARE) 


INPUT NO, OF PTS, 
1 Ss ) Ss 
INPUT LOWER LIMITS 


0. 2.0 = 
INPUT UPPER LIMITS 
2.0 2.0 


1041 POINTS RE 
2066 POINTS ac 
AVEMAGE= 


> Ls 
2.0 “2.0 -2.0 
2.0 2.0 2.0 
JECTED 
CEPTED 
7.€638 


MINIMUM ORAG CUNFIGURATJON 


co 
6.114940 0. 


KIND FS ABASE AMIN 
KIND oF S,ABASt + AMIN 
KIND oF S.ABASE SAMI 
KINDeFS.ABASE S AMIN 
KIND» FS,ABASE + AMIN 
KIND oF Se ABASE e AMIN 
KIND sFS.ABASE CAMIN 
KIND oF S.ABASE eAMIN 
KIND eFS.ABASE AMIN 
KINO of SeABASE SAMI 
KIND +FS+ABASE* AMIN 


KIND eFS.ABASE S AMIN 


VOLIIMF = RASELINE«M 


L 8 
-1.000000 

& 50.000 7465.94 
4 90.000 1003.> 
4 190,00 1438.0 
4 190,00 1641.8 
# 200,00 1665.3 
# 205,00 1706.5 
~~ 210,00 1747.9 
+ 215,00 1782.1 
« 220,00 1775.4 
& 225,00 2768.6 
4-06 272.37 
& 230,00 1761.8 


2. 


IN. CONF esNOSEoee 264920,7 


F 
000000 
705432 
100S.7 
1437.7 
1650.0 
1697.6 
1720.9 
1765.1 
1801.7 
1796.0 
1790.) 
256248 


1784.2 


266021.9 


4 
-1.000000 


26088.12 


at 
1.000000 


FIGURE 26, MINIMUM VALUE RESULTS WITH VOLUME CONSTRAINT 
(5x 5 LATIN SQUARE) 


INPUT ORDER OF FT = 8 THRU 14 
10R.]TPRINT 10 1 


VOLUME CONSTRAINT+ VOLUME 116480.8 GREATER THAN BASELINE 


ADDITIONAL TERMS ARE. « 2 
4105 4H 


INPUT DATA.... WAVE ORAG® 
1 ~1046020E-01 
& =. 9495500E-02 
7 .llee290e-01 
2  =«1105250€-01 
S  .1040130€-01 
8 .1142970E-01 
3... 1078680E-01 
6 .1263200€-01 
9 .1143090€-01 
0 .9793590E-02 
COEFFICIENTS OF wAVE ORAG EQUATION 
CH -10G0130E-01 ¢ 2) ~2700607E=03 (3) »S78S167E-03 » 3638333€-03 «7459000E-03 
6) +2598000E-03 « 7) 25391000€=03 ( 6) =21181679E-03 °2227321€-03 26171 788E-06 
STO DEV= .6¥863E-15 3 
FITTED WAVE OwAG 
eb s106K020E-01 ¢ 2) 21165250€-01 3) -10786R0E-01 ©9455500E-02 ~1060130€-01 
( 6) 1283200E-01 (¢ 7) 21166290E-01 8) 21142970€-01 211463090E-01 9793590E-02 


1 


ERROR oF FIT 
(1) *.2775558E-15 (2) -,6661538E-15 3) -22775558€-15 = 26661338E-15 -.1609823E-16 
(6) =.6106227E-15 ¢ 7) =, 304S761E=15 8) = ,6661338E-15 = .1665339E-15 = .4996006E-15 


POP e ee Hee Hee Hees eeHese Meet seretOttoes 


INPUT DATA. .e6 *VOLUME® 
1 =70402.30 
4 +-48181.60 
7 32094.00 
2 1205.500 
S 36291.10 
8 334646.10 
3 #7920.10 
6 1374659.1 
9  194575.0 
10 Oe 
GLEY 36291.10 48555. 30 71076.02 7140,667 (5) 14605,67 
( 6) 4258,>267 2179-207 8759.46) 8386.327 (10) -6849,6466 
STO DEV= .32d92€-08 
( 1) -79602.30 ( 12652500 «7920.10 -48181.60 (5) 36291,10 
(6) 1374659.1 ( 32094,00 83446.10 154575.0 (10) =,7275958E-09 


ERROR OF FIT 
(1) .2793968E-08 ( 2) OL - +2328306€-08 22095476E-08 (5S) -,1396984E-08 
(6) =,6519258E-08 ( 7) -,6924919E-09 -.3725290€-08 =25587935E-08 (10) -,727S958E-09 
SOO e EHH HOST HEHEHH ESSE SEHH eee etetHeose 
POSSESS eT HSS STHHeSH Hee tH esse tet esas 
INPUT NO. OF PTS. 
Zi 2k) 6 EI 
INPUT LOWER LIMITS 
-1.0 -1.0 -1.0 “1.0 
INPUT UPPER LIMITS 
1.0 wr) 1.0 1.0 
1792R2 POINTS REJECTED 
15199 POINTS ACCEPTED 
AVERAGE= .11728E-01 


MINIMUM ORAG CONFIGURATION 


co H A R1 R2 
«107408SE-01 26000000 - 29000000 -.8000000 «9000000 


116686.8 120002-2 


FIGURE 27, SURFACE FIT AND MINIMUM VALUE SEARCH RESULTS WITH 
VOLUME CONSTRAINT (3 x 3 LATIN SQUARE) 


orn opcadeen 


SECTION IV 


OPERATIONAL ASPECTS 


CORE AND TIME REQUIREMENTS 


Both the LATIN] and LATIN2 can be run on 60,000, words of memory. The 

computer time required for the search is dependent on the total number of search 
cases. As a guide it takes about 5-10 seconds for the LATIN1 program and 0,05- 
0.10 seconds for the Latin2 program to compute 1000 search cases. Additional 5 


seconds should be allowed for compilation. 
ERROR MESSAGES 


If the number of input data is less than the requested order of Latin Square fit 


(number of terms in the wave drag equation) the error message, 
*INSUF FICIENT NO, OF INPUT DATA* 

will be printed and the job will be terminated. 

VOLUME DATA PREPARATION 


If the volume constraint option is used, the two minimum value search programs 
require the input of configuration volumes as well as wave drag coefficients. For 

the LATIN] program, the input volume is that of the fuselage nose ahead of 

fuselage station 50, while for the LATIN2 program, the input volume is that of 

the complete configuration which comprises the fuselage, canopy and wings. Two 
auxiliary programs, which use the same body discription data for the 3DMoC program, 
can be used to calculate the input volume data. The listings of the two volume cal- 
culation programs, one for LATIN] and the other for LATIN2 programs, are pre- 
sented below. 


One input-data card is required for each of the volume calculation programs. The 


data required are 


KASE, Y, Y2 FORMAT (I5, 2F10, 2) 


where KASE = an input integer to identify the body configuration (e.g., cell number) 


Yl, Y2 = the fuselage stations between which the volume is to be calculated. 


The values of Y1 and Y2 used in these two programs are 


[rx rar] 
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3NNILNOD 


CCT) dxZ¢ (TD) KZOAS(TI9SS (TIL 8 (1)S* (1) 68 (1) 08 (1)d) SDINOD 1109 


2*tsI » 00 € 


(0°25*8°ST3S*2T) LYWeO 4 
(SIT*S°OT4) LVWYuO4 
3NNILNOSD 


B=RLW ((8°D39° 1) °YO* (6°03°1)) 4T 


€t=w (2°19°1) 4I 


(1)9S¢ (1) 29 (105901048 (100801) do 7 (10292) GV 3a 


N*Tt=y 2 00 

L=w 

N®&XVWA (0022) 034 
XVWA= IXYWA 

2 ONIM3SS (CIXVWAPL IEA) ST 


((IXVWASIOPA) PONV? (XVWAS3TP9A) ) ST 


0S*t=I 2 00 
Az(V1)SA 
VIK@AVILTSIO* TAZA 
TeVI= Vly 
AN*T=11 0S OG 
NX/ (TA*#2A) =AV1130 
ITeAN=NX 
ZA*STA*3ISUM (84S) 034 
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0022N 
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SECTION V 
LOGICAL STRUCTURE 


1, INTERDEPENDENCE OF SUBROUTINES 


The Calling-Called matrix for the program is shown in Figure 28. 


bo 
° 


LISTINGS 


Two complete listings of the Surface Fit and Minimum Search Program are 
given below, following Figure 28. The first listing is for the computer program 
developed in Phase I using a 6 x 6 Latin Square. The second listing is for the 


computer program developed in Phase II using a 4 x 4 Latin Square. 
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(a) LATIN1 Program 


(b) LATIN2 Program 


FIGURE 28, CALLING-CALLED MATRIX FOR SURFACE FIT 
AND MINIMUM SEARCH PROGRAMS 
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SECTION VI 


FUNCTIONAL DESCRIPTION OF SUBROUTINES 


XYGEN - Generates the Latin Square Coordinate Matrix. 
CROUT - Uses Crout method to solve simultaneous linear equations. 
SEARCH - Searches for the minimum value. 
VALUE - Uses the Latin Square Surface Fit Wave Drag Coefficients. 
XFORM - Transforms reduced variables to physical variables for forward 
fuselage. 
SIGMA2 - Performs integration using the second order Lagrange formula. 
CNSTN - Imposes constraints according to input parameter. 
BASEA - A second entry of CNSTN to transfer the value of the baseline 
cross-sectional area. 1 
VOLMD - A third entry of CNSTN to transfer the value of the fuselage 
volume, 
AREA - Computes the cross-sectional area of fuselage excluding the 
canopy. 
KONIC - Computes fuselage widths at given fuselage stations. 
KONICA - A second entry of KNOIC to compute the cross-sectional area 
‘ of the lower fuselage. 
TABLE - Compares the fuselage width distribution with given local 
C tabulation of width or with randome contour. 
‘ READ - Reads baseline body description data. 
4 
CONICF - Computes Y and Y' from given conic coefficients. 
LINEAR - Linear interpolation routine. 
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Appendix 


SAMPLE DECK SET-UP 


Samples of deck set-up are given below, The IVS and 3DMoC programs, 
because of their lengths, are stored on tapes in an updatable binary form. Permanent 
files c these programs are created for computation. The Surface Fit and Minimum 
Search Programs are read in the computer to create permanent files for computa- 
tion. The samples were selected to illustrate various calculations for wave drag 


reduction of an F-4 type baseline configuration. 
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